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The status of lattice calculations of heavy-light decay constants and of the B parameter 
Bg is reviewed. After describing the lattice approach to heavy quark systems, the main 
results are discussed, with special emphasis on the systematic errors in present lattice 
calculations. A detailed analysis of the continuum limit for decay constants is performed. 
The implications of lattice results on studies of CP violation in the Standard Model are 
discussed. 

1. Introduction 

Heavy quark systems have attracted a lot of interest in the past decade. Theo- 
retical and experimental investigations of weak decays of heavy quark systems offer 
the main sources of information on the phenomenon of CP violation in the Standard 
Model. Our understanding of CP violation is directly related to the knowledge of the 
elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, which describes the 
pattern of quark mixing in flavour-changing charged current interactions. Within 
the framework of the Standard Model the CKM matrix is unitary, and, provided 
that one can determine its elements with high enough precision, any deviation from 
unitarity would be; interpreted as a hint of "new physics" . The theoretical treat- 
ment of weak decay amplitudes involving heavy quarks is, however, hampered by 
large uncertainties due to strong interactions. Since quarks are confined within 
hadrons, the exchange of soft gluons between them makes weak decay matrix ele- 
ments intractable in perturbation theory. Therefore, other theoretical tools have to 
be applied in order to deal with the intrinsic non-perturbative nature of weak decays 
of hadronic systems, e.g. QCD sum rules ^ or lattice gauge theories Another 
formalism which can be applied specifically to heavy quark systems is provided by 
the Heavy Quark Symmetry and the Heavy Quark Effective Theory (HQET) 

Lattice gauge theories offer the possibility of a systematic non-perturbative 
treatment of strong interaction effects in weak decay amplitudes. Indeed, Monte 
Carlo simulations of lattice QCD have made important contributions to our un- 
derstanding of hadronic physics by providing non-perturbative, model-independent 
estimates for a number of quantities, which are relevant for the spectrum and weak 
decays of hadronic systems in general. An important ingredient in non-perturbative 
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studies of weak decays is the operator product expansion . In general, a weak am- 
plitude, which usually involves the product of two hadronic currents, can be written 
as a sum over local operators, viz. 

^(f-i)=^C„(M)(f|6„(/x)|i). (1) 

n 

Here, fj. is the renormalisation scale, and the Wilson coefficients C„ (/i) , which con- 
tain the short-distance physics, are independent of the initial and final states and 
can be computed in perturbation theory for /i 3> Aqcd- The matrix element 
(f|0„(/x)|i), on the other hand, describes the long-distance physics, contains strong 
interaction contributions and must be calculated non-perturbatively, e.g. in a lattice 
simulation. 

The role of lattice simulations in the analysis of weak matrix elements is two- 
fold. First, they serve to test QCD, as lattice data can be compared to experiment 
or to results obtained using HQET or QCD sum rules. Second, they can make 

predictions for yet unmeasured quantities like the decay constant Jb- In fact, fs 
(or rather the combination /bVBb, where .Bb is the B parameter relevant for 
B^-B^ mixing) is the principal unknown quantity in the analysis of CP violation 
in the heavy quark sector. If lattice QCD is to have a serious impact on testing the 
consistency of the Standard Model, then it is clear that the control over systematic 
errors in the lattice approach is of paramount importance. 

There exists by now a vast amount of literature on lattice results for heavy- 
light decay constants, which will be discussed in detail in this review, pointing out 
the differences between various calculations and assessing their inherent systematic 
errors. Another aim of this study is to combine as many lattice results as possible 
in a systematic fashion, and to provide the currently best estimate for fBy/Bs 
computed on the lattice. The present status of lattice calculations of B meson 
decay constants and B parameters can be summarised as follows 

fs = 172 i^^MeV, /b,//b, = 1.14 ±0.08, (2) 

fB\fK = 195 t^°MeV, BB = l.Stl, (3) 

where -Bb is the renormalisation group invariant B parameter. In the following 
sections we will review lattice data from several groups and describe the derivation 
of the global results above. 

An overview of the development of this area of lattice gauge theories can be 
gained by consulting various review talks in the proceedings of recent lattice confer- 
ences Qtiigp reviews and pedagogical introductions can be found in refs. i9~26_ 

We begin this review in section 2 with a brief introduction to lattice QCD for 
the non-specialist and describe the lattice approach to heavy quark systems. In 
sections we discuss lattice results for heavy-light decay constants. A detailed, if 
somewhat technical analysis of the continuum limit is presented in subsection 3.4, 
which may be skipped if the reader is not so interested in the technical details. 
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Section 4 contains the discussion of lattice results for the B parameter Bb- The 
implications of our findings on the study of CP violation is presented in section 5. 
Finally, section 6 contains some concluding remarks. Two appendices are added, 
which present details of the matching between lattice matrix elements and their 
continuum counterparts, as well as a summary of lattice data in the light hadron 
sector. 

2. Lattice Approach to Heavy Quark Systems 

In this section we present the basic concepts of lattice QCD before various formu- 
lations of heavy quarks on the lattice are discussed in some detail. Furthermore, the 
relation between weak matrix elements on the lattice and their continuum counter- 
parts is discussed, and at the end of this section we summarise the main systematic 
errors which affect lattice results for heavy quark physics. 

2.1. Lattice QCD and numerical simulations 

Since its original foundation by Wilson in 1974^'^, lattice QCD has developed 
into a mature area of research in elementary particle physics. In general, the lattice 
approach now forms an important part in the study of a large class of quantum field 
theories, ranging from scalar field theory to quantum gravity. 

One important aspect in the study of lattice field theories is that they are suitable 
for numerical methods such as Monte Carlo simulations. Indeed, the biggest impact 
lattice QCD has had on phenomenology in the past decade has originated from large- 
scale numerical simulations. The strength of this approach is based on the fact that 

the lattice is currently the only known non-pcrturbative regulator of a quantum 
field theory which is systematically improvable. 

In this subsection we will briefly discuss the basic features of the formulation and 
the implementation of lattice QCD in Monte Carlo simulations. For more detailed 
information, the reader is referred to several good textbooks 

In order to be easily accessible to computer simulations, lattice QCD is usually 
formulated in euclidean space-time. For a finite system, the lattice size is ■ T, 
where is the spatial volume, and T is the extension of the lattice in the time 
direction. The lattice sites are separated by the lattice spacing a. It is easy to 
show that the lattice spacing acts as an UV cut-off. Present computers can handle 
lattices of typically 16* to 48* lattice sites. Thus, to be able to accommodate a 
hadron, the physical length of the lattice should not be smaller than about 1.5 fm. 
This implies that values of the lattice spacing in physical units must lie in the range 
of a = 0.1 0.05 fm. 

On the lattice, non-abelian gauge fields are represented by so-called link vari- 
ables, which connect neighbouring sites. These link variables, which carry a Lorentz 

index fi, arc elements of the gauge group (here: SU(3)). Given a link variable U^{x) 
emanating from site x in direction /U = 1, . . . , 4, one can define a vector field Afj^ix) 
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as an element of the Lie algebra of the gauge group SU(3) via 

giaA^(x) = g gu(3). (4) 

The simplest lattice version of the action of a Yang-Mills theory is the Wilson 
plaquette action^ 

Sg[U] = /3E E (l - ^ReTrP^.(x)), /3 = 6/gl (5) 

where Pn„{x) = Ufj,{x)Uy{x + jl)Ujj^{x + i')Ul{x) is the product of links around 
an elementary square of the lattice, the "plaquette", and is the bare gauge 
coupling. Here, fi denotes a vector in direction /i of length a, such that a; + /t is 
the neighbouring site of x in direction jx. Using the definition of the field Afj,(x) in 
eq. (4), one can show that for small lattice spacings, the term in Sg\U] reduces to 
the familiar expression — ^TrF^j,. 

Quark and antiquark fields, ip{x), i^ix) are associated with the lattice sites. 
Using the link variable J7^(a;), a possible definition of the covariant derivative is 
given by 

D^i,{x) = }-{U^{x)i,{x + fi)-i,{x)] (6) 
V^(a;)5^ = -^{^{x + il)Ul{x)-^{x)}. (7) 
In euclidean space-time, the Dirac matrices can be defined to satisfy 

{7^,74 = 2V- (8) 

Now we are in a position to write down a latticised version of the QCD Lagrangian. 

It turns out, however, that the naive lattice transcription of the fermionic part 
suffers from the notorious fermion doubling problem: in the case of a free Dirac 
particle it simply means that, in four dimensions of space-time, the propagator 
contains 16 poles, resulting in a 16-fold degeneracy of the spectrum. Following 
Wilson's proposal^, the degeneracy can be lifted by adding a counterterm (i.e. 
the Wilson term) to the action which has the effect of pushing the masses of the 
unwanted doublers to the cut-off scale. However, the Wilson term breaks chiral 
symmetry explicitly, which has wide-ranging consequences as we shall sec below. 

Another method for the removal of doubler states is the use of "staggered" or 
Kogut-Susskind fermions^°. In this approach, the individual spin components of a 
Dirac spinor are spread over several lattice sites within a hypcrcube. Thereby the 
number of degenerate states is reduced from 16 to 4, which are then interpreted as 
different flavours. It should be emphasised, however, that the spin/flavour assign- 
ment is only valid in the continuum limit. This is different for Wilson fermions, 
where spin/flavour assignments on the lattice are exactly as in the continuum. This 
makes the identification of local operators much easier, and therefore the Wilson 
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formulation appears to be better suited for the study of weak matrix elements of 
local composite fields. However, staggered fermions leave a U(l) x U(l) subgroup 
of chiral symmetry invariant, and therefore the use of staggered fermions may, after 
all, be advantageous for studying quantities such as the B parameter Bk, whose 
chiral behaviour is a central issue. As most results for heavy quark physics have 
been obtained using Wilson quarks, we shall concentrate on this approach in the 
remainder of this review. 

The Wilson action for lattice QCD reads 



where 



S[U, V] = So [U] + [U, V, V] , (9) 
f * 1 

X ^ fl — 1 

+^(x + /i)(r + 7^)f/^(x)V(x)] +V^(a;)V(a;)|, (10) 
= a''J2^ix)M{x,y)i>iy), (11) 

and in the last line we have introduced the Wilson-Dirac operator A4. The Wilson 
term is the piece proportional to r in the above expression. The hopping parameter 
K is related to the bare mass mo via 

i^=7i (12) 

2amo + 8r ^ ^ 

and the Wilson parameter r is usually set to one. The Wilson action is thus conve- 
niently parametrised in terms of bare parameters (/3, k) instead of the bare gauge 
coupling and quark mass {go, mo). The parametrisation of the fermionic part of 
the action in eq. (10) implies that the quark and antiquark fields have undergone a 
rescaling according to 

ip{x) V2Kip{x), i){x) ip{x) \P1k. (13) 

This normalisation is referred to as the relativistic norm of quark fields. 

The addition of the Wilson term proportional to r to the naive lattice action 
leads to an additive renormalisation of the quark mass. This implies that there 
exists a critical value of the hopping parameter, Kcritj at which the quark mass 
vanishes and chiral symmetry is restored. The subtracted quark mass is defined by 

2 ^ AC K^crit ^ 

and in the free theory the critical value of k occurs at 

1 , , 

Kcrit = o' 50 = 0. (15) 
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The role of a ^ as an UV cut-off can also be seen from the fact that at small 
coupling 

5o'~l/lna, P = 6/9l (16) 

and therefore the continuum limit, a — > 0, is formally reached as ,3 — > oo. Note 
that typical values of the lattice spacing in physical units are ~ 1.5-4.0 GeV in 

current simulations, which corresponds to /3 ~ 5.7-6.4. 

By means of the lattice discretisation procedure, one has suitably altered the 
UV (short-distance) behaviour in order to obtain a non-perturbatively regularised 
theory, whilst preserving the long-distance physics. The continuum result for, say, 
hadron masses computed on the lattice is then obtained in the limit a ^ 0. At this 
point it is worth noting that the Wilson fermion action differs from the classical 
continuum action by terms of order a. In contrast to this, the leading discretisation 
errors (lattice artefacts) for the pure gauge action arc only O(a^), and there arc 
arguments that this is also true for staggered fermions"^^. The predictive power 
of lattice QCD depends crucially on the degree to which lattice artefacts can be 
controlled. This is of particular importance for heavy quarks, where the effects due 
to finite lattice spacing are large, as we shall sec below. 

The lattice action £'[[/, in eq. (9) can now be used to define a partition 
function 

Z = J D[U]D[i^]D[il;]e-^^^''^''f'\ (17) 

where, on a finite lattice 

/ D[U] =11 [ dU,{x), (18) 

and dUiJt{x) is the invariant group measure. The vacuum expectation value (VEV) 
of an observable O is defined as 

(0) = 1 y"D[C/]£)[V^]£)[V']C»e-^°[^l-^^[^''^'^l 

= ^ j D[U\0 detM[U]e-'^''^^\ (19) 

where in the last equation we have introduced the determinant of the Wilson- 
Dirac operator by integrating out the fermion fields. The aim of any Monte Carlo 
simulation is to evaluate (O) stochastically through an average O of individual 
measurements 0{Ui}, computed on a set of gauge configurations {Ui}, i = 1,. . . ,Nc 

- 1 

{0)c.O=—^Om, (20) 

i 

where a sequence of gauge configuration has been generated in a Markov process 
using det A^[f] e~'^°[^l as probability measure. This implies that Monte Carlo 
results have a statistical error that decreases proportionally to 1/ \/iVc. In the limit 
of infinite statistics, the relation between (O) and O becomes exact. 
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If the determinant det A^[/7] is to be included in the probabihty measure, then, 
roughly speaking, it needs to be evaluated in every update step in the Markov chain. 
Since det M [U] is highly non-local this is prohibitively costly on present computers. 
Various algorithms have therefore been constructed which should allow for a more 
efficient evaluation but still the computational effort is very large. 

It has therefore been proposed to set det A^[[/] = 1 in the generation of gauge 
configurations, which corresponds to neglecting effects due to quark loops. This 
procedure defines the so-called quenched approximation which is still most 
widely used in current simulations. Consequently, the overwhelming part of the 
material covered in this review was obtained using the quenched approximation. 
As far as hadronic quantities are concerned, the quenched approximation amounts 
to computing observables in QCD in the presence of a Yang- Mills background field. 

We end this subsection with a few more technical remarks. Hadronic observables 
are extracted from mesonic or baryonic two- or three-point functions, which are 
constructed from quark propagators Sq{x,y). The quark propagators themselves 
are obtained as the inverse of the Wilson-Dirac operator M. by solving 

M{x,y)Sg{y,z) = 5^,. (21) 

Hence, on every gauge configuration {Ui\, i = 1, . . . ,Nc, one needs to invert the 
matrix that couples the quark fields in the lattice action. Therefore, the use of an 
efficient inversion algorithm is of great importance. The most widely used inversion 
algorithms'^ include the Conjugate Gradient, Minimal Residual and, more recently, 
the Stabilised Biconjugate Gradient method''*. 

As an example of a mesonic two-point function constructed from quark prop- 
agators, we now discuss the euclidean correlation function of the axial current 
= '07^75?/', which is used to extract properties of pseudoscalar mesons. For 
euclidean times t = X4 > 0, the correlation function is defined by 

C{t;p) ^ e-^^~-^"(0|A4(a-,i)4(0)|0), (22) 

X 

which, after performing the Wick contractions, can be written in terms of quark 
propagators Sg{y,x) as 

Cit;p) = e-'P^-(TV{745,(0,ar)745t(0,a;)} ). (23) 

Here we have used that Sq{x,Q) = 755*^(0, a;)75. Inserting a complete set of inter- 
mediate states in eq. (22), one obtains the spectral decomposition of the correlation 
function, viz. 

Cr(,;a = j:!(i^M!e-«' (24) 

For large t, the lightest state in the spectral decomposition dominates. For p = 0, 
the asymptotic form of C{t;p) is given by 

C(,Ol-^°K^^#Ke---, (25) 
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where Mp is the mass of the pseudoscalar meson in the ground state. 

In order to obtain reliable estimates of the mass Mp and the matrix element 
(01^4(0) I P) in a lattice calculation, one has to be able to follow the signal for 
C{t; 0) up to large t, which requires good statistics. Furthermore, the lattice has to 
be long enough for the lightest state to be observable, and therefore one normally 
chooses T = 2L. For heavy quark systems, the issue of reaching the asymptotic 
behaviour of C{t; 0) is a particular problem; here, the correlation function falls off 
rapidly and may have dived into the statistical noise before the ground state can 
be observed. Therefore, one has to choose operators in the evaluation of C'{U 0) 
which have a good overlap onto the ground state. Using point-like, local opera- 
tors such as V(a^)7^75V'(a;) above, results in a very poor signal. This can be intu- 
itively understood from the fact that the wave function of the particle one wants 
to study is an extended object and not point-like. Therefore, when solving for the 
quark propagators, eq. (21), one applies so-called smearing techniques, which are 
designed to generate more extended sources, having a better overlap onto the de- 
sired state Since the singal-to-noise ratio for hadron masses and 
decay constants can be enormously improved using "smeared" sources and/or sinks 
in the propagator calculation, smearing techniques are now an indispensable tool in 
lattice simulations of QCD. 

In a simulation, dimensionful quantities like hadron masses or pseudoscalar de- 
cay constants are obtained in lattice units. In order to convert these lattice estimates 
into physical units, one needs to set the lattice scale. This is usually done by com- 
paring a low-energy hadronic quantity computed on the lattice to its physical value. 
For instance, using the mass of the p meson one obtains [GeV] via 

where aMp denotes the lattice result for the p mass. Other quantities used to set 
the scale include the pion decay constant f^, the nucleon mass, the string tension 
^/a■ or the hadronic scale ro discussed by Sommer''^. 

One complication that arises in the computation of aMp (and, indeed, of all 
quantities that involve u and d quarks), is that one cannot simulate directly at the 
u and d quark masses. This comes from the fact that hadronic quantities involving 
very light quarks have large correlation lengths, which would lead to severe finite- 
size effects on cm-rently accessible lattice sizes. Therefore, one typically uses quark 
masses in the region of that of the strange quark and extrapolates the results to mu,d, 
or to the chiral limit. Finite-size effects are discussed in more detail in subsection 2.4. 



2.2. The b quark on the lattice 

As already mentioned in subsection 2.1, current simulations are carried out at 
typical values of the inverse lattice spacing in physical units of ~ 1.5-4 GeV. 

For heavy quarks, this implies that one needs to worry about the effects of finite 
cut-off already when one wants to study charm physics, since TOcharm is only about a 
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factor of two smaller than . Furthermore, it is evident that the b quark cannot be 
studied directly, since its Compton wavelength is smaller than the lattice spacing. 
Hence, one concludes that lattice artefacts are large for heavy quarks, and that 
special care is required to control them. 

In the following we describe several methods which arc employed to circumvent 
the problem that the b quark cannot be studied directly or are designed to alleviate 
the problem of large discretisation errors. 

2.2.1. 0{a) improvement 

In this approach one seeks to reduce lattice artefacts by using so-called im- 
proved actions and operators in order to cancel the leading discretisation error in 
on-shell quantities. The concept of perturbative improvement was first outlined by 

Symanzik""' and developed further in rcfs.^^^^^. For the Wilson action in lattice 
QCD, Sheikholeslami and Wohlert have shown that the 0{a) discretisation error in 
the action can be cancelled by adding a local counterterm^^ 

Ssw[U,^,i^] = Sg[U] + 5f [/7,?, V] + csw^a' ^ ii{x)a^,F^4x)i;{x), (27) 

where the improvement coefficient Cgw depends on the gauge coupling. In order 
to obtain 0(a) improved matrix elements of composite operators, these composite 
fields have to be improved, too**^'^°. For instance, the improved axial current reads 

A^ix) = A^,{x) + CAad^P{x), (28) 

where ^^(a;) = tp{x)jf_i'^5'ip{x), P{x) — tl;{x)-f5^p{x), and ca denotes another im- 
provement coefficient. Recently, the ALPHA Collaboration has proposed an 
0{a) improved action for which the improvement coefficients Cgw and ca were de- 
termined non-perturbatively for (3 > 6.0. This action has so far not been used in 
simulations of heavy quark systems. 

2.2.2. Non-relativistic normalisation 

In a second approach, it has been suggested to suitably adapt the Wilson 
action, such that the Wilson propagator does not deviate from the continuum be- 
haviour even for quark masses am 1 in lattice units (i.e. for quark masses above 
the cut-ofl[). It has been argued that this can be achieved by a modified rescaling 
of lattice quark fields (c.f. eq. (13)) according to 

^(x) ^ \/2^e"'"p/V(a;), (29) 

where the "pole mass" amp of the Wilson propagator is given by 

amp = ln(l + am), (30) 

and am is defined in eq. (14). The factor \/2Ke'""''/^ is designed to interpolate 
smoothly between the relativistic and non-relativistic regimes. 
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Consciquently, in order to cancel the effects of large quark masses in studies of 
matrix elements of composite operators, the normalisation of quark fields is changed 
according to the above scale factor. This non-relativistic normalisation of quark 
fields according to eq. (29) is referred to as the Kronfeld-Lepage-Mackenzie (KLM) 
norm. 

Higher-order corrections for large quark masses can also be considered. As is 
argued in ^^'^^'^6, this can be motivated by the observation that for quark masses 
am > 1, the mass amp and the "kinetic mass" am2 appearing in the non-relativistic 
dispersion relation, are no longer equal 

aE(p) = amp + h . • . , m2 ^ mp, (31) 

2am2 

and the relevant quark mass is the kinetic mass TO2. In the free theory, the kinetic 
mass am2 is obtained from the free Wilson propagator in the non-relativistic limit 

eamp sinh(amp) 

am2 = . , , \ — (32) 

smh(amp) -|- 1 

Both 0{a) improvement and the KLM norm arc used for quark masses in the 
region of that of the charm quark. Clearly, residual lattice artefacts remain and 
must ultimately be extrapolated away (which, however, can be performed much 
more reliably if they are small). The results obtained around mcharm must also 
be extrapolated to the mass of the b quark, and clearly one needs to control this 
extrapolation in order to obtain meaningful results. In what follows, we shall refer 
to this approach, where relativistic heavy quarks are used in conjunction with a pre- 
scription to reduce lattice artefacts (either 0{a) improvement or the KLM- norm), 
as the "conventional method" . 

2.2.3. Static approximation 

Here the b quark is treated as infinitely heavy Prom an expansion of the heavy 
quark propagagtor in the inverse heavy quark mass, l/mg, one obtains at leading 
order ^ 

SQ{x,t;0,0) = |e(t)e-™«*i±^ + e(-t)e'"«*i^| S{x)Vs{t,0), (33) 

where V(^{t, 0) is the product of links from ((), t) to the origin, for example for < > 0, 

pQ{t, 0) = ul{6, t - l)Ul{6, t-2)--- C/i(0, 0). (34) 

This form of the propagator is also obtained using a Lagrangian for heavy quarks 
Q{x) which simply reads 

£static^Qt(a,)£,^Q(a.)^ (35) 



where the covariant derivative is defined in eq. (6), and Q(x) is a two-component 
spinor describing the heavy quark. 
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One expects corrections of order Aqcd / fng to the results in the static approxi- 
mation, which can potentially be large. With a few exceptions, such as the B*-B 
mass splitting, the computation of higher order corrections in l/mq to the static 
limit is complicated due to the presence of power divergencies^^. Nevertheless, the 
static approximation is a valuable tool in lattice studies of heavy quark systems. 
It plays the crucial role of guiding the extrapolation of results obtained with the 
conventional method to the mass of the b quark by providing direct information at 
infinite quark mass. 

2.2.4. Non-relativisUc QCD (NRQCD) 

A non-relativistic formulation of heavy quark systems can be obtained by ex- 
cluding relativistic momenta through the introduction of a finite cut-off Auv ~ mq 
such that 

p rnqv <^mQ, (36) 

where v is the 4- velocity of the heavy quark. The loss of relativistic states through 
this requirement can be compensated for by adding new local interactions order 
by order in p/Auv- In this way one obtains a cut-off theory as an expansion of 
the original action in p/Ajjy ~ v. The non-relativistic QCD Lagrangian is then 
obtained by applying a Foldy-Wouthuyscn transformation, which separates quark 
and antiquark fields, thereby generating an expansion of the QCD Lagrangian in 
l/mg. For instance, at order l/mg the NRQCD Lagrangian reads 

^NRQCD ^ Qtfj^^ _ ^)q- Q^^-^Q. (37) 

\ 2mQ J 2mQ 

This approach is thus based on the a priori discretisation of the non-relativistic 
formulation of the Wilson action. Although the above expression contains the term 
that defines £>^*3,tic gq (35), one should realise that there are important differences 
between NRQCD and the static approximation: the theory defined by the NRQCD 
Lagrangian is non-renormalisable, since a finite cut-off has to be kept 

Auv~a~^^TOQ- (38) 

In the language of the lattice this means that the formal continuum limit, a — > 0, 

does not exist, in contrast to the case of the static approximation. Hence, for 
NRQCD to work in a lattice simulation, one needs to calculate at fairly large values 
of a (i.e. at small /3). Lattice artefacts (i.e. cut-off effects) have to be reduced by 
including higher orders in l/mg. For the study of matrix elements of composite 
fields using NRQCD to nth order in l/niQ, it is important to take into account 
corrections of ©(I/ttiq) in the operators as well. 

Prom the above discussion of various methods to formulate heavy quarks on the 
lattice, it is obvious that none of them is entirely satisfactory. They are subject to 
rather different systematic effects and thus provide complementary information on 
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heavy quark systems. The fuU picture wiU therefore only emerge when results from 
all methods are compared. 

2.3. Weak matrix elements on the lattice and their continuum counter- 
parts 

It has already been noted in subsection 2.1 that the Wilson action in eq. (10) 
breaks chiral symmetry explicitly. Therefore, even in the massless theory, vector 
and axial vector currents are not conserved, since the UV behaviour of the theory 
has been changed as a consequence of the regularisation procedure. In other words, 
the regularisation procedure conflicts with the naive conservation of the currents. 
The resulting short-distance corrections between the vector and axial currents in 
lattice and continuum theories can be absorbed into normalisation factors Zy and 
Za, respectively. 

The chiral Ward identities for the axial current normally ensure that the ax- 
ial current does not get renormalised. However, since chiral symmetry is broken 
explicitly by the Wilson term, the Ward identity no longer applies, as it is vio- 
lated by terms of order a. Nevertheless, the condition that the correctly normalised 
lattice axial current satisfies the Ward identities can be used to derive its normali- 
sation Since the normalisation factors incorporate short-distance effects, they 
can in principle be calculated in lattice perturbation theory. 

Using the unimproved Wilson action the one-loop perturbative results for the 
vector and axial current normalisation constants read''^'''^'''^ 

ZA = l-0.n3gl + O{g^), = 1 - 0.174 .g^ + 0(g4)^ (39) 

whereas for the 0{a) improved Sheikholeslami-Wohlert action* eq. (27) one ob- 
tains 

Za^ I ^0.018 g^ + 0{g^), Zy = I - O.lOO g^„ + 0{gt). (40) 

These expressions are written in terms of the bare gauge coupling g^ , which is given 
by gQ =6/(3. However, it has been argued in^'' that the bare gauge coupling is a bad 
expansion parameter due to the appearance of large giuonic tadpole contributions 
in the relation between the link variable [/^ (.x) and the continuum gauge field (x) 
(see eq. (4)). The poor convergence property of lattice perturbation theory can be 
improved by absorbing these tadpole contributions into suitable redefinitions of the 
bare parameters of the theory. The simplest definition of an improved expansion 
parameter is the so-called "boosted" coupling g"^ introduced by Parisi 

g'^i, (41) 

where Uq is taken to be the measured average value of the plaquette (^ReTrP) 
{uo = (|ReTrP)^/* is a gauge invariant estimate of the average link). Replacing 

*Note that here we quote the expressions for the non-local currents which have been "rotated" 
according to the prescription given in 
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the bare coupling by in the one-loop expressions for Za and Zy defines the 
so-called "boosted" perturbation theory. 

In the fermionic sector, one can absorb tadpole contributions into a redefinition 
of the hopping parameter «; in a similar fashion according to 

K = KUo- (42) 

It is then expected that the critical value of k is much closer to its tree-level value 
of 1/8, such that uq — l/8Kcrit) c.f. eq. (15). This is the basic concept of tadpole 
or mean field improvement^^. By using non-perturbative input for uq such as 
the measured link or. alternatively, Kcrit, one can improve the UV behaviour of 
the lattice theory, which should manifest itself in a better convergence of lattice 
perturbation theory. 

Tadpole improved expressions for the normalisation factors are obtained by com- 
bining their perturbation expansions with those for uq as well as the measured value 
of the latter For instance, the tadpole improved one- loop expression for the axial 
current normalisation reads 

Za = uo[i + [Z^2^ -ui'^]g^), (43) 

where Z^^ and Uq'^ are the one-loop coefficients in the expansions of Za and uq, 

respectively. The numerical value of u^^"' depends on whether one identifies uq with 
the average link or the critical hopping parameter, thus considering the perturbative 
expansions of either (iReTrP)V4 or SKcrit- 

In conjunction with the KLM-norm defined in eq. (29), tadpole improvement 
also plays an important role in reducing lattic;e artefacts due to large quark masses. 
As an example we consider the case of the lattice axial current V'i7ai75V'2 of quark 
fields tpi, tjj2 with hopping parameters Ki and K2, respectively. The normalisation 
factor of the axial current is then supplemented by an extra factor of 

exp{a(mp,i -I- mp,2)/2}, (44) 

which arises from the modified rescaling of quark fields in eq. (29). Inserting the 
definition of the pole mass and applying tadpole improvement, the normalisation 
factor becomes 



,a(n^P,+™p,)/22^ = J(l + :^-^)(l + ^-:rJ-)ZA 

yV 2ki 2Kcrit^\ 2k;2 2Kcrit^ 

= \/(SH(SH^- 

where in the last step we have used Kcrit = ^critWo — 1/8, a relation that becomes 
exact if uq is defined by 1 /S/^crit rather than the average link. 

Apart from the coupling constant defined in eq. (41) other definitions of a mean 
field improved coupling can be used. Lepage and Mackenzie proposed a coupling 
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av defined through 

ay(3.41a-^)(l - (1.19 + 0.017n/)ay) = In(iReTrP). (46) 

If av is to be used for the evaluation of a tadpole improved normalisation factor, 
then, according to ref. it has to be computed at a scale q*, which denotes the 
mean momentum flow relevant for a given matrix element. Typical values of q* lie 
in the range 1 < aq* < tt. 

Given the potential uncertainty in the perturbative expressions from higher or- 
ders, it is desirable to perform non-perturbative determinations of the normalisation 
factors. For Za and Zy this can be done by imposing the chiral Ward identities 
as a normalisation condition A more general approach, which can be ap- 

plied to a large class of operators ®^ , imposes the normalisation condition between 
quark states. Numerical values have so far been obtained for the axial current for 
Csw = 0^^, in the static approximation for Cgw = 1 and also for the AS* = 2 four- 
fermion operator relevant for the Kaon i? parameter Bk^'^'^^- Systematic errors 
in non-perturbative determinations of normalisation factors are, of course, present 
and need to be controlled. Recently, the ALPHA Collaboration has determined 
Za and Zy non-perturbatively in the whole range < < 1, with total errors at 
the 1% level. 

Apart from the problem of the normalisation of lattice operators, another conse- 
quence of explicit chiral symmetry breaking is the possibility of mixing of operators 
with a definite chirality. For instance, the AF = 2 four-fermion operator defined 

by 

Ol = (^7^(1 - 75)V') (^7^1 - 75)V') (47) 
can mix with its right-handed counterpart 

Or = (V^7^(l + 75)V') (V^7/.(l + 75)V') • (48) 

Furthermore, operators can mix with higher dimension operators, which is also 
a consequence of the breaking of Lorentz invariance by formulating QCD on a 

hypcrcubic lattice. 

Hence, in general the relation between the matrix element of an operator O in 
the continuum and in the Wilson formulation is given by 

{f\6\ir-' = ^Z„{f\&:''\i) +0{a), (49) 

a 

where a labels the operators of the same dimension that can mix with each other, 
and Za are the short-distance normalisation factors. 

2.4. Summary of systematic errors 

The main systematic effects that affect lattice results for heavy quark physics 
can be broadly divided into lattice artefacts, finite-size effects, quenching errors and 
normalisation of lattice operators. 
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In any simulation in physical volumes, the following inequalities must be satisfied 

a «: M-i «: L, (50) 

where M is a typical hadronic mass, and L denotes the spatial length of the lattice. 

The left part of the inequality constrains the possible values of heavy quark masses 
that can be studied safely, i.e. without suffering from large discretisation errors. 
On the other hand, very light quark masses will give rise to finite-size effects, as is 
evident from the right part of the above inequality. Hence, a major limitation of 
current lattice simulations is that very different scales need to be incorporated on 
a single lattice. 

We now list and discuss the main systematic errors which affect lattice results 
for heavy quark systems: 

2.4.1. Discretisation errors 

Discretisation errors, i.e. the effects of the finiteness of the lattice spacing, are 
present in all simulations of lattice QCD, but are especially important for heavy 
quark systems as we have seen in subsection 2.2. We have already discussed how 
0{a) improvement can be used to reduce these effects, and there are attempts to ex- 
tend the improvement programme beyond 0{a) ^^'^^ or to even construct "perfect" 
lattice actions that ideally are completely free of discretisation errors Once 
these cut-off effects are reduced, either by employing the Symanzik improvement 
programme or by absorbing effects of large quark masses into the normalisation 
of quark fields (KLM-norm), residual discretisation errors have to be extrapolated 
away by taking the limit a ^ 0. 

It should be emphasised that this procedure does not solve the problem that the 
b quark cannot be studied directly in current simulations using relativistic quarks. 

2.4.2. Finite-size effects 

Lattice estimates of hadron masses and matrix elements are distorted due to the 
finiteness of the lattice volume. An analytic calculation by Liischer showed that 
the mass shift in finite volume decreases exponentially with L 

6M{L) = M{L) - M(oo) ~ e-^l^°. (51) 

Numerical studies, however, revealed that for spatial lengths of Z/;$2.0fm the ob- 
served mass shift deviates from the exponential decrease and is very well described 
by a power law behaviour. Using a simple model, the effect can be ascribed to 
the "sqeezing" of hadronic wave functions for small volumes and for values 
oi L^2 fm, one recovers the exponential behaviour of eq. (51). 

In practice one obtains estimates of hadronic observables in infinite volume by 
repeating the simulations at different values of L and extrapolating to the infinite 
volume limit. 

Finite-size effects also prevent light quark propagators from being calculated at 
realistic values of the masses of u and d quarks, whose bound states have large 
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correlation lengths. Therefore, light quark masses are usually taken around the 
mass of the strange quark, and hadronic observables are then extrapolated to the 
chiral limit, using the quark mass dependence deduced from chiral perturbation 
theory as a guide for the extrapolation. 

2.4.3. Quenching 

As most lattice results for heavy quark systems are still obtained in the quenched 
approximation, attempts must be made to quantify the effects due to neglecting 
quark loops, in order to provide meaningful information for phenomenology. In 
the past few years, several collaborations have reported results for heavy quark 
physics using dynamical quarks, but much more precise simulations have yet to be 
performed. 

An indirect manifestation of quenching effects is the uncertainty in the; lattice 
scale [GcV]. This uncertainty arises because different quantities used to set the 
scale give different results, which is understood from the fact that quark loops make 
different contributions to the various quantities used to compute [GeV]. 

2.4.4. Normalisation and mixing of lattice operators 

Uncertainties in the normalisation factors relating lattice matrix elements to 
their continuum counterparts come from the fact that these factors are, in most 
cases, only known perturbatively to leading order in g^. In a few cases, e.g. for the 
heavy-light axial current in the static approximation, this uncertainty may be of the 
order of 10%. Non-perturbative estimates of these constants are therefore highly 
desirable. For vector and axial vector currents the chiral Ward identities can easily 
be employed as normalisation conditions in a non-perurbative determination. For 
more complicated operators or scale-dependent renormalisations, a more general 
condition needs to be applied ®^ , which, however, may be less accurate numerically. 

In order to check the consistency of different methods, one needs to make contact 
between pcrturbativc and non-pcrturbativc; estimates of normalisation constants. If 
gQ is used as an expansion parameter in the evaluation of the perturbative expres- 
sions, this is only possible at large values of /?, where non-perturbative methods 
normally fail. Using a "boosted" or tadpole-improved coupling g^ instead, one may 
try to make a comparison in a region where current simulations arc being performed. 
In some cases it was found that "boosted" perturbation theory yields values close 
to the non-perturbative determinations, but the issue remains inconclusive as long 
as small values of the couplings cannot be investigated numerically. A systematic 
non-perturbative analysis for < gig < 1 has been performed in the framework of 
the Schrodinger functional 

3. Leptonic Decays of Heavy Mesons 

In this section we review lattice results for heavy-light decay constants such as 
fo, /ds, /b and Jb,- After introducing the basic definitions, we present and discuss 
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the results obtained using different formulations of heavy quarks on the lattice. 
Subsection 3.4 contains a detailed discussion of the approach of the continuum limit 
for unimproved Wilson fermions. Readers who are not interested in the details 
of the continuum extrapolation will find the main results listed in subsection 3.5, 

where also a comparison to other theoretical methods is presented. Throughout 
this review we use a convention for which f^^ = 131 MeV. 

3.1. Basic definitions 

In the continuum, the pseudoscalar decay constant fp is defined by 

(0|^40)|P(p)) = zp^fp, (52) 
whereas the vector decay constant 1/fv is usually defined by 

{0\V,mV) = e^-^. (53) 
Jv 

Hence, in order to extract the decay constants of heavy-light mesons, one needs to 
compute the matrix elements of the axial and vector currents between a mesonic 
state and the vacuum. 

On a lattice, this is achieved by analysing various mesonic two-point functions 
of the form 

C£1W^E(0|'^?(^)4''(0)|0), (54) 

X 

where Ji and J2 are interpolating operators which can annihilate or create the 

heavy-light pseudoscalar or vector meson under study. The superscripts Q, R de- 
note whether a local (L) or smeared (5) interpolating operator is used. 

One possibility to extract the heavy- light pseudoscalar decay constant is to con- 
sider the ratio 



cm) ~ j:Ao\p'ix,t)p^'{om 

|^-(M.(r/.-<,), 

where we have used the expressions for the asymptotic behaviour for both C\p{t) 
and Cp^(t) on a finite lattice with time extension T and periodic boundary con- 
ditions. Using the smeared pseudoscalar density P^{x) to create the pseudoscalar 
meson in Cap and Cpp results in a much better signal for these correlation func- 
tions. This is of particular importance in the static approximation discussed below, 
where it is notoriously difficult to obtain a reliable signal, and where calculations 
using purely local operators are known to fail^^'^^'®^'^^^. 

In order to determine (OIA4 (0)|P) which contains the decay constant, one also 
needs to know the pseudoscalar mass Mp and the matrix element (0|P'^(0)|P), 
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which are extracted from separate fits to 

C|f (t) mP'mP)\'' g-MpT/2 ^^^^ (Mp(r/2 - t)) . (56) 
2 Adp 

From (0|^4 (0)|P) one can then extract fp via 

{0\Ai{0)\P)=Mpfp/ZA, (57) 

where Za is the normahsation constant of the axial current discussed in subsec- 
tion 2.3, which is usually evaluated in ("boosted" or tadpole-improved) perturba- 
tion theory. The actual value of Za depends on whether an improved action has 
been used (i.e. on the choice of Cgw), on the details of the implementation of tad- 
pole improvement, and also on whether the relativistic, non-relativistic or static 
formulation of heavy quarks has been used. 

In the two-spinor formalism of the static approximation (see eq. (35)), there is 
no distinction between A4{x) and P{x), such that the pscudoscalar decay constant 
is extracted from combinations of the following correlation functions 

C^^it) = J2(^\Ai{x,t)Al'mO)*>^° {A^fe-^^ (58) 
C^^{t) = ^(0|Af(a;,f)Af (0)|0) '^5° A'^A^e-^*, (59) 

X 

where £ is the unphysical difference between the meson mass and the mass of the 
heavy quark. The decay constant /^^* is related to the matrix element A^ via 

y^tat ^ j^L ^2/ Mb Zf""* . (60) 

The one-loop expressions for Z^'** were calculated in refs. 59,60,6i^ g^j^^j later extended 
to the 0{a) improved case^^'^^. Here we do not list the expressions but refer 
the reader to appendix A, where we have listed the normalisation factors for some 
operators of interest. 

A similar procedure as the one described above can be applied to determine the 
vector decay constant 1/ fv- 



3.2. Decay constants and heavy quark symmetry 

Lattice results for decay constants using the conventional or non-relativistic 
formulations can be used to test predictions of the Heavy Quark Effective Theory. 
It is well known that HQET predicts a scaling law for heavy-light decay constants 
in the limit of infinite quark mass, mg oo 

^-fp-'-^a-y^o^ (61) 

where Mp ^ My ^ M ^ mg. Hence, the heavy quark flavour symmetry implies 
that fp^Mp behaves like a constant (up to logarithmic corrections) 

' const a-'^"^\ (62) 
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where (3o = 11 — and n/ is the number of active quark flavours. The heavy 
quark spin symmetry predicts the vector and pseudoscalar decay constants to be 
degenerate. Defining U{M) as the ratio of matrix elements of the axial and vector 
currents one expects 

U{M) = - 1, mg ^ ^. (63) 

In order to study the mass dependence of fp\/Mp, one divides out the scaling factor 
Q^'^^^° and defines the quantity 

$(Mp) ^ /pv/m; f^^^ V''" • (64) 



as{MB) 

In figure 1 we compare some results for <i?(Mp) obtained using the static approxima- 
tion with those from the conventional formulation One observes 
rather high vahies in the static approximation, whereas the results using rclativis- 
tic heavy quarks are much lower. The figure also illustrates that improvement 
is a crucial ingredient for detecting deviations from the scaling law: using either 
the Sheikholeslami-Wohlert (SW) action in eq. (27) or employing the KLM-norm, it 
seems possible to interpolate between the static approximation and the conventional 
approach. Using the standard relativistic norm in eq. (13) instead leads first to a 
flattening and then to a decrease of $(Mp) as the heavy quark mass is increased, as 
is shown in figure 1 using the results of the PSI-CERN-Wuppertal (PCW) group 
Hence, the infinite mass limit of $(Mp) would be in complete disagreement with the 
results obtained directly at infinite h quark mass. Figure 1 also illustrates our earlier 
remark that different formulations of heavy quarks have to be used simultaneously 
in order to reveal the full picture. Another remarkable observation is that results 
from simulations with different systematics agree very well within errors (at simi- 
lar values of the lattice spacing), which underlines the consistency of the different 
methods. 

From the slope of $(Mp) one can infer the size of 1/Mp corrections to the 
scaling law, eq.(62). They amount to about 15% at the mass of the B meson and 
about 40% at the mass of the D. Thus it appears that there are large corrections 
to the predictions of the heavy quark flavour symmetry for decay constants. Recent 
studies using NRQCD i04-i08 j^aye reported even larger corrections at Mb- 

A test of the heavy quark spin symmetry can be performed by computing the 
quantity U{M) and studying its mass dependence. In the heavy quark limit, U{M) 
differs from one by short-distance corrections 

UiM) - ^ = (l + + 0(1/M)) . (65) 

The UKQCD collaboration performed this analysis using the 0{a) improved SW 
action; the short-distance corrections were divided out, and the resulting data were 
fltted to either a linear or quadratic function of the inverse, spin-averaged mass 
1/M. 
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Fig. 1. The quantity $(Mp) plotted against the inverse pseudoscalar mass. Diamonds represent 
the data from ref. obtained at ,8 = 6.3, squares denote the results from UKQCD loi.ioo^ using 
the SW action at /3 = 6.2. The data from refs. 98'102 at /3 = 6.26 are represented by plus signs 

(KLM-norm) and crosses (rclativistic norm). The data obtained in the static approximation are 
slightly shifted for clarity. Also the positions of the B and D meson masses are indicated. 



Figure 2 shows that U{M) indeed is consistent with one in the heavy-quark 
hmit. This result demonstrates not only the manifestation of the heavy-quark spin 
symmetry in the infinite mass limit, it also provides support for the parametrisation 
of the non-scaling behaviour of fp\/Mp as a power series in 1/Mp. 

Despite these encouraging results, heavy quark scaling laws for decay constant 
require further investigation. In particular, the size of 1/M corrections has to be 
studied as one approaches the continuum limit, in order to detect to what extent 
the slope in fp^Mp is influenced by lattice artefacts. Also, it is a priori not clear 
whether the interpolation of fp^/Mp between the results obtained in the static limit 
and the convc;ntional approach is feasible at non-zero values of the lattice spacing. 
Since discretisation errors are potentially very different in the two approaches, a 
meaningful interpolation may only be possible after extrapolation to the continuum 
limit. We will return to this point in subsect. 3.4. 

3.3. Results for pseudoscalar decay constants 

We now compile and discuss the results for the pseudoscalar decay constants 
Jd, fs, fos and Jb^ from various authors and using different formulations of heavy 
quarks on the lattice. 

3.3.1. The conventional approach 
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Fig. 2. Data for the quantity [/(M) from the UKQCD collaboration . The solid and dashed 
lines represent the linear and quadratic fits to the data. Radiative corrections to U{M) have been 
divided out. 

We start with the decay constants fo and fo^ , which have only been calculated 
using the conventional (relativistic) formulation. The determination of fo proceeds 
by interpolating the lattice data for $(Mp) (c.f. eq. (64)) to the mass of the D me- 
son, as either a linear or a quadratic function of 1/Mp. In table 1 we list the results, 
also indicating the value of /3, the normalisation of the quark fields and whether an 
improved action was used (i.e. Cgw > 0). The first error on the decay constant is 
the statistical error, whereas the others are an estimate of systematic errors. The 
MILC collaboration ^^^^^ ^ ^^[j.^ 

error which is an estimate of quenching 
effects, and which is obtained by comparing the results in the quenched approxima- 
tion at ~ 2GcV to a simulation with n/ = 2 flavours of dynamical quarks at 
approximately the same value of a~^. 

Table 2 shows the corresponding results for Jb and Jb, ■ Furthermore, the SU(3)- 
flavour breaking ratios fosl fo and fBjfs are listed in table 3. The extraction of 
Jb using the conventional method involves the extrapolation of $(Mp) in 1/Mp 
to the B system, which is an additional source of uncertainty. In order to have 
better control over this extrapolation, some authors prefer to interpolate between 
the values obtained in the static approximation and those from the conventional 
method As discussed in subsection 3.2, at a fixed value of (3 this 
appears only possible if either an improved action or the KLM norm is employed. 
Using instead the relativistic norm with Csw = results in a flattening of $(Mp) as 
Mp is increased, which is interpreted as a lattice artefact. This is consistent with the 
observation that early calculations using the conventional method obtained rather 



22 Leptonic decays of heavy quarks on the lattice 



Table 1. Results for fo and fos using the conventional formulation of heavy quarks. All values 
were obtained in the quenched approximation, except for those from the HEMCGC collaboration 
which were computed for Uf = 2. Infinite values of f3 symbolise results obtained after extrapolation 
to the continuum limit. Results by MILC and JLQCD are preliminary. 
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low values for fn and fs i26.i25^ a,nd underlines the importance of addressing the 
effects of large quark masses in the conventional approach. 

In addition to using the KLM norm in conjunction vi^ith the otherwise unim- 
proved Wilson action (Cg^ = 0), some collaborations ^^'ii°'iii'ii2'ii5 also include 
higher-order corrections for large quark mass, as explained in subsect.2.2. The 
pseudoscalar meson mass Mp has been shifted accordingly, i.e. 

Mp ^ M'p = Mp + (m2 - mp), (66) 

where mp and m2 are defined in eqs. (30) and (32), respectively, and also tadpole 
improvement has been applied tomp andm2. The mass dependence of $(Mp) is 
then studied in terms of the shifted mass Mp, which is the tree-level estimate of 
the kinetic mass appearing in the non-relativistic dispersion relation. 

Tables 1 and 2 show that the results among different collaborations are broadly 
consistent, provided that the effects of large quark masses have been treated, either 
by employing an improved action, or by using the KLM norm. The error analyses by 
the different groups reveal that the statistical errors amount to 5-10%, so that the 
uncertainty in pseudoscalar decay constants is now dominated by systematic effects. 
In view of the large overall errors and the different treatment of the systcmatics, 
it would be premature to conclude on the basis of tables 1 and 2 that discretisation 
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Table 2. Results for fs and /s^ using the conventional formulation of heavy quarks. All values 
were obtained in the quenched approximation, except for those from the HEMCGC collaboration. 
Results by MILC and JLQCD arc preliminary. 
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Table 3. Results for the SU(3)-flavour breaking ratios fo,/ fo and fB,/fB using the conventional 
formulation of heavy quarks. All values were obtained in the quenched approximation. Results by 
MILC are preliminary. 
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errors are under control. This requires a detailed analysis of the continuum limit, 
which we shall attempt in the next subsection. 

The ratios JdJ/d and /bJ/b listed in table 3 vary by about 10% among differ- 
ent groups, which is slightly larger than the typical statistical error. This is inspite 
of the fact that some of the systematic errors in these ratios (e.g. the renormalisa- 
tion factor Za) are expected to cancel, such that they are determined much more 
reliably. 

Here we wish to make a few comments about the estimation of systematic er- 
rors. Most collaborations quote a systematic error coming from the uncertainty in 
the lattice scale. Apart from the statistical error in the quantity that is used to set 
the scale, this may also include the effects of choosing different quantities to esti- 
mate a~i [GeV]. The MILC Collaboration ii2'ii3 include this uncertainty in their 
estimation of quenching errors. Furthermore, MILC perform a detailed analysis of 
systematic errors after the extrapolation to the continuum limit. Their estimate 
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yields a total uncertainty in (b duo to using large quark masses of about 10% of the 
result. Further sizeable systematic effects are ascribed to the extrapolations in the 
lattice spacing and the light quark mass. By comparing results for /b for n/ = 
and n/ = 2 at similar values of [GeV], MILC conclude that /s may be larger 
by around 10% in the unqucnchcd theory. 

Bhattacharya and Gupta have analysed different prescriptions in the defini- 
tion of the tadpole improved normalisation factor Za- Using either uq = (|Re Tr P) 
or uo = 1/8/tcrit in cq. (45) leads to an uncertainty of 14 MeV in fjj, which is consid- 
ered a conservative estimate by the authors. They also ascribe an error of around 
5% to the uncertainty in fixing the charm mass, using either the kinetic mass M'p or 
the conventional "pole" mass Mp for the heavy-light pseudoscalar meson. However, 
since only a single mass value for the heavy quark is used in the authors are not 
able to perform a detailed analysis of the mass behaviour, and thus the problem of 
using different definitions of heavy masses cannot be fully addressed in their study. 

In fact, AUton et al. have found that the mass dependence of $(Mp) is only 
marginally changed by considering the "kinetic" mass instead of the "pole" mass. 
They quote their best estimates from their run at /? = 6.2, Cgw = 1, using results at 
/3 = 6.0 as well as data with Cgw = to estimate systematic errors. Furthermore, 
AUton et al. emphasise that the data for <I>(Mp) using Csw = 1 are indistinguishable 
from those for Csw = computed with the KLM prescription, in good agreement 
with ref. and this work (c.f. figure 1). 

JLQCD investigate mass effects over a range of quark masses and lattice spacings 
using the unimproved Wilson action At fixed values of (3 they analyse the 
effects of using different field normalisations (relativistic or KLM) and definitions of 
the pseudoscalar mass ("pole" or "kinetic"). Whilst different prescriptions produce 
a fair amount of variation in the value of the decay constant at non-zero lattice spac- 
ing, JLQCD's preliminary findings indicate that these deviations are decreasing 
as a ^ 0. Similar observations were made earlier by the PCW collaboration 

3.3.2. The static approximation 

Now we turn to the discussion of results for Jb obtained in the static approx- 
imation. In table 4 we list the results from several simulations. In general, f§^^ 
turns out much larger than the results from the conventional approach. Especially 
some of the earlier simulations produced very large values which may 

partly be due to the fact that the axial current normalisation factor in the static 
approximation, Z^^*", was evaluated in bare perturbation theory, thus producing a 
much larger value. 

In order to avoid discussing many different systematic effects in fg^^ when com- 
paring different simulations, we have also listed the matrix element in lattice 
units in table 4 (see eq. (60)). By comparing results at fixed /?, one observes con- 
sistency in among many different collaborations over a wide range of /3 values. 
This is remarkable as even the use of an improved action for the light quark docs 
not seem to have a discernible effect on A^ (see also the discussion in^°°). Thus, 
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whilst the numerical values for obtained using either Cgw = or Cgw = 1 are 
totally consistent at fixed /3, the one-loop estimates of the renormalisation factor 
^stat g^j.g lapggj. 10-15% for Cgw = 1 as compared to Cgw = 0. Therefore, as long 
as there is no non-perturbative determination of Zj***, which may clarify the issue, 
lattice determinations of fg^^ will always differ in the unimproved and improved 
theories at non-zero values of a. It is therefore of great importance to study the 
approach to the continuum limit of both cases. 

As the signal-to-noise ratio in simulations of the static approximation is notori- 
ously bad ^''''''1'^^, it is important to use an efficient smearing technique in order to 
extract a reliable signal. The most advanced calculations in the static approxima- 
tion 130,128,100 variational approach in which a matrix correlator is constructed 
by computing cross correlations from different smearing functions. One can then 
diagonalisc the matrix correlator and project onto the ground state, thereby elimi- 
nating approximately the contamination of the signal from higher excited states 



Table 4. Results for fg, fg^/fB and the matrix element in lattice units, a^^^^A^, obtained in the 
static approximation. All data are from quenched simulations. The results in ^'^'^ are preliminary. 



Collab. 


/3 


/|tat [MeV] 


IbJIb 


a3/2A^ 




UKQCD 


6.2 


01^1^ -t-18 -t-28 
^DD _20 _27 


1.16 


0.112 tl 


1 


APE 


6.4 


235(9) 


1.13(2) 


0.075(3) 


1 




6.4 


209(14) 


1.12(7) 


0.076(5) 







6.2 


221(12) 


1.12(7) 


0.109(6) 


1 




6.1 


190(10) 


1.13(5) 


0.135(7) 







6.0 


258(9) 


1.19(3) 


0.201(7) 


1 


FNAL128 


"oo" 


188(23)(15) t'^(14) 


1.22(4)(2) 









6.3 


225(17)(14) 


1.17(3)(1) 


0.099(8) 







6.1 


215(21)(14) 


1.23(3)(2) 


0.135(13) 







5.9 


241(13)(13) 


1.21(2)(1) 


0.250(14) 







5.7 


271(13)(20) 


1.18(3)(1) 


0.564(28) 





APE 129 


6.2 


290(15) (45) 


1.11(3) 


0.111(6) 


1 


UKQCD 101 


6.0 


nop, + 8 +67 
ZOD _iQ _42 


1.13 ti 


0.211 


1 




6.3 


235(20)(21) 


1.11(2)(2) 


0.092(6) 





Ken 130,131 


6.0 


224 tr 


1.22(1) 


0.184(7) 





SH132 


6.0 


297(36) 1^ 




0.206(25) 





APE 133 


6.0 


370(40) 


1.19(5) 


0.22(2) 


1 




6.0 


350(40) (30) 


1.14(4) 


0.23(2) 





PCW98 


"oo" 


230(22)(26) 


1.16(5) 







PCW118 


6.0 


366(22)(55) 


1.10(9) 







HS 134 


6.0 


386(15) 









ELC 135 


6.0 


310(25)(50) 




0.22(2) 






A possible interpretation of the very low value for observed by the Ken- 
tucky group 130, who use a large basis of different smearing functions, is that there 
may be some residual contamination from higher excited states in the results of 
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Tabic 5. Results for /s using NRQCD. The results by the SGO Collaboration 10B,106 ^g^e obtained 
on dynamical configuration with Uf = 2. 



Collab. 






/b [MeV] 


Ib, I Ib 




SGO 
SGO 


6.0 
5.6 




183(32)(28)(16) 
126-166 


1.17(7) 
1.24(4)(4) 



2 


Ken 


6.0 


0.110 tt 




1.15 tl 





gjjl32 


6.0 


0.156(8) 


171(22) til 








refs. ^^^'^^^ °°. This is not implausible, since typical values for the binding energy 
extracted from the exponential fall-off of the correlation function in eq. (58) are 
around £ ~ 0.7 in lattice units at /? = 6.0. This is still rather large and thus great 
care is required in isolating the ground state. Therefore, despite using smeared 
sources, and in view of large overall statistical fluctuations, the results at the lower 
end of the range in /3 could still be affected by higher excited states if no variational 
approach is applied. The situation improves above (3 > 6.0, where the binding 
energies are lower such that not so much can be gained by employing variational 
smearing. At /3 = 6.2, UKQCD have explicitly verified that the contamination 
from higher states was negligible, and that their results with and without varia- 
tional smearing were consistent 

The dependence of on the number of dynamical quark flavours n/, has 
been studied in refs. The results by the Rome 2 group indicate that 

at = 1.1 GcV the value of /I'*** increases by 15 20% for nj = 3. Further 
calculations are clearly needed to confirm this sizeable flavour dependence. 

3.3.3. Non-relativistic QCD 

So far there are relatively few estimates of fs using NRQCD. For all but one 
study an estimate for the axial current normalisation constant in NRQCD was 
not available, such that most previous studies had resorted to using Z^'***. The 
major results are listed in table 5. Quantities in which the current normalisation 
cancels, such as Jb, / /b are consistent with the results from the static approximation 
and the conventional approach. 

A major advantage of NRQCD in studies of the scaling law for fpy/Mp is that 
it is possible to study the corrections in the inverse meson mass 1/M starting from 
the static approximation as the limiting case. Indeed, very large 1/M and 1/M^ 
corrections to the scaling law eq. (62) have been observed in NRQCD 
However, as the NRQCD hamiltonian in those studies was obtained after trunca- 
tion at order I/tuq in the heavy quark mass, one may question the observation 
of large 1/M^ corrections. More recent studies have therefore included all I/toq 
terms in the NRQCD action and operators. Preliminary results by Ali Khan and 
Bhattacharya indicate that the individual corrections of 1/M^ to the NRQCD 
hamiltonian amotmt to about 20% of the 1/M corrections, although their effect on 
the slope of fpy/Mp is marginal, as the various 1/toq corrections to the hamiltonian 
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and operators come in with opposite signs. 

All of these findings must, however, be re-evahiated in the light of a recent study 
in which Za computed in NRQCD at one-loop order, has for the first time been 
applied. This has a dramatic effect on the 1/M corrections to the scaling law, which 
turn out to be much smaller when Za is taken into account. Furthermore, the scale 
q* at which the mean field improved coupling a{q*) in the one- loop expression for 
Za is evaluated (c.f. subsection 2.3) has a big influence on the mass behaviour of 
fp\/Mp. Choosing q* = a^^, /pV Mp is almost constant in 1/M, whereas for 
q* = Tr/a the slope of jp\JMp in 1/M is consistent with the findings using the 
conventional approach Performing a careful estimation of systematic errors, 

the authors of ref. quote 

/s = 174(28) (26) (16) MeV, = 1/a, (67) 

/b = 183(32)(28)(16)MeV, g* = 7r/a, (68) 

where the first error combines the statistical and the fitting error, the second comes 
from the uncertainty in the lattice scale, and the third is an estimate of uncertainties 
due to higher-order contributions to Za and neglected higher orders in 1/M. 

There are indications that consistency in the mass behaviour of ]p\fMp from 
the region of charm to the static limit is observed using different formalisms of 
treating heavy quarks on the lattice. It should, however, be emphasised that it 
is crucial to investigate the continuum limit of the mass dependence of fp\/Mp, 
which we will attempt in the next subsection. 

3.4. The continuum limit 

We are now in a position to study the continuum limit of heavy-light decay 
constants using results discussed in the previous subsection. Although some of the 
more recent simulations 98ii02, 128,111, 112, 115 j^^ve already performed an extrapolation 
to the continuum limit, we want to treat as many different results as possible on 
an equal footing and perform a thorough analysis of systematic errors in the final 
results. 

As is apparent from tables 1-4, most results have been obtained for Wilson 
fermions with Cgw = 0. At present it is therefore not possible to perform the 
extrapolation a — > of the data with Cgw = 1, which would be an important 
consistency check. Furthermore, since the data of the MILC m-H^ g^j^j JLQCD 
collaborations are still preliminary and are likely to be updated, we leave them out 
in the following analysis, despite their high statistics and large range in /3. 

Since the systematic effects in the results listed in tables 1-4 are quite different, 
e.g. the choice of lattice scale, the details of the quark field normalisation and 
in the evaluation of Za, a straightforward extrapolation of all listed results would 
be unreliable. Therefore, in order to study the mass dependence of fp\fMp and 
obtain estimates for decay constants in the continuum, we adopt a similar strategy 
as described in^^'^°^. The various steps of the procedure are as follows: 
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(a) At a fixed value of P one forms the dimcnsionlcss ratio '^{Mp) / ^/a using 
lattice data for $(Mp) defined in eq. (64) and the string tension 

(b) The ratio $(Mp)/v^^^^ is then interpolated to common values of Mp/^. 

(c) At every fixed value of Mp/^/a one extrapolates ^{Mp) / ^^^'^ to the con- 
tinuum limit as a function of a^/a. 

(d) Now one sets the scale using extrapolated lattice data for, e.g. the pion decay 
constant, Jt^J \fa\a=Q of the p meson mass Mpj \fa\a=Q- The continuum value 
of $(Afp)/Y^' is interpolated to the desired value of the heavy-light meson 
mass, and finally one converts back to the decay constant in physical units 
using the definition of 3>(Mp). 

This procedure implies that in step (d) /b is obtained through an interpolation 
between the result in the static approximation and those using relativistic heavy 
quarks. Therefore, a separate extrapolation of the data for /^*** is required (i.e 
step(c) for Vo'/Afp = 0). 

The difference between this procedure and the extrapolations performed in 
refs. 128,112,115 jg ^jjg^^ ^j^g latter the conversion of heavy-light decay constants 
into physical units precedes the extrapolation to a = 0. This is feasible if the scale 
is measured as well at every value of (i considered. The universal treatment of a 
number of different simulations, however, forces one to use the extrapolated lattice 
scales, since the values of /3 at which the quantity to set the scale is measured do 
not necessarily coincide with the /3 values used in the computation of heavy-light 
decay constants. 

Lattice data for the string tension a^fu used in the following extrapolations are 
listed in appendix B. In order to perform a controlled extrapolation to the continuum 
limit of a number of hadronic quantities, it has been suggested'*^ that the hadronic 
scale ro is to be preferred over the string tension -^/a. However, despite the merits 
of choosing ro to monitor the approach of the continuum limit, it has so far only 
been determined in a fairly limited range of /? (see refs. i95- 198,202 -j^ Therefore we 
stick to the string tension for our purposes. 

Appendix B also describes the details of the continuum extrapolation of quanti- 
ties in the light quark sector, such as /^r, Mp and the IP-IS splitting in charmonium, 
Aip_is. The pion decay constant j-^ is widely used to set the scale, because both 
/tt and the heavy-light decay constant /p are subjected to similar systematic errors 
- at least in the conventional approach - such as quenching errors and uncertainties 
in Za- It can therefore be expected that some of these effects cancel if the con- 
version is performed using /^r. In rcf. ^""^ it was argued that the IP-IS splitting in 
charmonium is largely insensitive to lattice artefacts, which also makes it a reliable 
quantity to set the scale. The mass of the p meson, on the other hand, shows a 
strong dependence on the lattice spacing for unimproved Wilson fermions (see e.g. 
refs. such that it is not clear a priori whether or not one controls its extrapo- 
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lation to the continuum limit well enough. We include the mass of the p meson 
nevertheless in this study. 

In this review we choose to quote our best estimates using to set the scale. 
Systematic errors are estimated by varying 

• the details in the evaluations of Za and Z^^* 

• the field normalisation of propagating heavy quarks (KLM or relativistic) 

• the quantity that sets the scale 

• the minimum value (3^i^ which is used in the extrapolation in step (c) above 
(i.e. the maximum value of a [fm]) 

3.4.1. Extrapolation of the data obtained in the static approximation 

We start our analysis by considering extrapolations of /g^* for Cgw = 0, using 

the tadpole improved value for Z^'*^* (see appendix A) . In figTire 3 we have plotted 



the ratio /bV-^bP^'^^/v^ as a function of the string tension in lattice units. 




Taking the data from refs. it can be seen from the figure 



that the simultaneous extrapolation of all data is not feasible and results in a large 
value of x^/dof. We have therefore chosen to divide the data sample into two sets 
that each can be extrapolated with acceptable x^/dof. Combining the data from 



refs. 99,127,128,130,132 j^to "sct 1" and those from refs. 98'i32,i33,i35 jj^^o "set 2", we 



note that the values in set 2 arc in general higher than those in set 1. As discussed 
in the previous subsection, a possible explanation for this may be the presence of 
residual contamination from excited states, because no variational smearing was 
used in set 2. However, since it is impossible at this stage to further investigate the 
apparent discrepancies between sets 1 and 2, we will take set 1 for our best estimate 
for fg^*', whilst quoting the extrapolated value from set 2 as a systematic error. 

We now discuss the dependence on the minimal f3 value, /3min, used in the 
extrapolation. Figure 3 shows that the extrapolated result for /b-s/MbI'**'^*/^/^'^^^ 
from set 1 changes if the point at (3 = 5.7 is included or not. This may signal that 
higher corrections in the lattice spacing may be important if one includes data for 
which a can be as large as 0.17 fm. This observation is consistent with the findings of 
ref. ^'^^ , where it was concluded that the expected scaling behaviour of the quantity 
is observed down to values of /? ~ 5.9 — 6.0. However, since one obtains a 
perfectly good fit of set 1 if the result at (3 = 5.7 is included, we again quote the 
difference between results obtained for /3min = 5.7 and /Jmin = 5.9 as a systematic 
error. 

Our best estimate is computed for /3min = 5.9 and we obtain 





= 1.81 ± 0.16 (stat) 



+0.31 
-0.20 



(extr) ± 0.07(norm) 



(69) 



■3/2 
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Fig. 3. The ratio Jb \/Mb\ '''^'"^ / y^^^^ for Csw = plotted versus the string tension. Open squares 
denote the results from refs. 12S, 127, 99, 130, 132, 133, 98, 135 ("ggj wlicrcas open triangles represent 
refs. 98il33,l36,l32 ("get 2"). The solid line represents the linear extrapolation of the data in set 1 
for /3jnin = 5.9, the dotted line is the extrapolation of set 1 for /3min = 5.7, and the dashed curve 
denotes the extrapolation of the data in set 2. The tadpole improved normalisation factor of the 
axial current, Z^** was used. 

where the first error is statistical, the second arises from the difference between 
setsl and 2 as well as using /Smin = 5.7. The third error is due to the variations 
from using different prescriptions in the numerical evaluation of Z^'^^. 

Using fjr/\/^\a=o = 0.285(13) from the extrapolation of the data of the GFll 
collaboration ^"^^ (see appendix B) and /,r = 131 MeV to convert into physical units, 
we find 

fs'"'' = 245 ± 27 (stat) tfg (syst) MeV, (70) 

where we have combined different systematic errors in quadrature, including those 
from using Mp and Aip_is to set the scale. Given the large errors, this result is 
not incompatible with the fairly low result quoted by the FNAL group and that 
of the PCW collaboration^* (see table 4). It should be noted, however, that the 
result in eq. (70) is a consequence of omitting the results at low /3 (i.e. (3 = 5.7) 
in the extrapolation. If those data are included one obtains the lower value of 
'^t = 218 ± 20 (stat) MeV. 

3.4.2. Extrapolation of data using relativistic heavy quarks 

Now we proceed to discuss the extrapolation of heavy- light decay constants using 
the conventional approach. This will finally enable us to study the scaling law for 
fp\/Mp, eq. (62) in the continuum limit. 
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In order to be able to interpolate fp^/Mp to common values of Mp whilst 
comparing the use of the kinetic versus the pole mass, we need to resort to data 
whose mass dependence has been published (as a function of the hopping parameter 
of the heavy quark). Therefore we only use refs. 99>i02,i09,ii6,ii9^ leaving the data by 
MILC j^j^j JLQCD 114,115 -(-Q included in the future. Furthermore wc chose 

not to use the data of ref. As already noted in^"^, it is questionable whether the 
asymptotic behaviour of the relevant correlation functions has been reached in^''^, 
since no smearing has been applied in the computation of the relevant correlation 
functions. Failure to isolate the ground state leads to an overestimate of the matrix 
element which is used to extract fp. In fact, the quantity $(Mp) from ref. is 
very different compared to refs. when all data are scaled appropriately 

and treated equally. Also some of the very early calculations 121,125.126 apply 
smearing techniques to improve the isolation of the ground state, and consequently 

we do not include them in our analysis. 

0/0 

The data for fp\fMp arc used to compute ^(Mp)/y/a , which is then interpo- 
lated to fixed values of y/a/Mp between 0.12 and 0.29, corresponding to heavy-light 
meson masses of 1.4-3.5 GeV. In ref. only a single value of the heavy quark mass 
was considered, which corresponds to an inverse pole mass of \falMp ~ 0.283 or 
to ^/a/M'p ~ 0.260 if the definition of the kinetic mass, eq. (66), is used. Only at 
these values of the inverse pole and kinetic masses, respectively could the results of 
ref. be included in the extrapolation to the continuum limit. 

In order to test the stability of the extrapolations at fixed s/a/Mp, we compare 
four procedures for which different evaluations of Za and different concepts to deal 
with heavy mass effects are compared: 

(1) KLM norm, kinetic mass, tadpole improved Za 

(2) KLM norm, pole mass, tadpole improved Za 

(3) relativistic norm, pole mass, tadpole improved Za 

(4) relativistic norm, pole mass, Za in boosted perturbation theory 

Here, the kinetic mass is always evaluated using eq. (66), in contrast to some results 
in refs. ii^'^o^ -^^iiere the kinetic mass was extracted from the measured dispersion 
relation. It should be mentioned at this point that using the KLM norm together 
with the pole mass may, after all, be an ill-defined concept as large mass corrections 
are only addressed at leading order, thus resulting in an incomplete treatment of 
these effects. In figure 4 we plot the extrapolation of $(Mp) to o = for all four 
schemes listed above. A remarkable feature is that for y/ajMp > 0.21, the results 
in the continuum limit are rather consistent, including the data computed using the 
KLM norm in conjunction with the pole mass. Since at fixed /3 the data computed 
using the relativistic norm show the wrong qualitative behaviour (see figure 1), 
this is quite a surprising result. Thus, at this stage it appears that the details 
of the procedure employed to deal with large mass effects and the normalisation 
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Fig. 4. The quantity ^(Mp) / ^/a' '~ plotted as a function of the string tension in lattice units. 
Open squares represent the data computed procedure (1) as described in the text. Full triangles 
correspond to procedure (2), and the crosses and open circles represent the data using proce- 
dures (3) and (4), respectively. The solid lines denote the linear extrapolations of '■P{Mp)/\/a' ' 
to the continuum limit, using only the data for which /3 > 6.0. The dotted lines are the continuation 
of the fits to the point at /? = 5.74. 



factor Za produce a variation of only 3-5% in the continuum results for $(Mp) for 
meson masses Mp ^ 2.2 GeV. These observations are consistent with the findings 
in refs. However, the spread of results at a = using different procedures 

increases for larger masses, as can be seen in figure 5. This is an indication that 
eventually lattice artefacts make the continuum extrapolation harder to control if 
the quark mass is chosen very large. The slope of $(Mp)/V^^/^ as a function 
of Qy/a increases for larger heavy quark mass, which is the expected behaviour as 
the influence of lattice artefacts becomes greater. As figure 4 shows, the results 
computed using the KLM norm and the kinetic mass show the smallest slope in the 
lattice spacing a compared to all other procedures employed. Unlike in the case of 
the static approximation, the extrapolations of ^{Mp) / ^/a for relativistic heavy 
quarks are fairly insensitive to the value of /?min- This can be seen from figure 4 where 
the data for ^{Mp)/^/a'' at (3 = 5.74 are not incompatible with the extrapolations 
performed for /3min = 6.0. However, below we will still use /3i„in = 6.0 for our best 
estimates as a safeguard against lattice effects in the final results, especially since 
the lattice spacing at /3 = 5.74 is only ~ 1.1 GeV. 

3.4.3. Interpolation to the physical meson masses 

3/2 

We are now in a position to perform step (d) and interpolate ^{Mp)/^/a 
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Fig. 5. The quantity $(Mp)/y^^''^ in the continuum hmit as a function of \/7t/Mp. Open 
squares denote the data computed using procedure (1), whereas the crosses represent the points 
obtained using procedure (3). Both sets were obtained using f^min = 5.7 in the extrapolation. For 
^{Mp)/i/a^^^ in the static approximation, the circles denote the extrapolated values for set 1 
using /3niin = 5.7 (open symbol) and /^min = 5.9 (full symbol). The value obtained using set 2 is 
denoted by the triangle. 



to the physical meson masses. Setting the scale by fTr/\/^\a=o = 0.285(13), we 
have interpolated our results for ^{Mp)/^/^^^^ at y/a/Mp = 0.21, 0.24 and 0.26 
to y/a/Mo = 0.245(11). Our final result for fo thus is 

fo = 191 ± 19 (stat) t 5(extr) tjolscale) MeV, (71) 

where the first error is statistical, the second is an estimate of systematic effects 
by varying /Jmin and the treatment of large mass effects. The third error is due to 
choosing Mp or Aip_is to set the scale. The latter effect is predominantly due to 
the known mismatch of fjr/Mp computed in the quenched approximation and the 
experimental result, which is partly ascribed to quenching. 

In the case that the light quark has been interpolated to the mass of the strange 
quark, steps (a)-(d) of the extrapolation procedure are entirely analogous, and for 
we obtain 

fo, = 206 ± 17 (stat) t 7(extr) (scale) MeV. (72) 

By interpolation of the data for ^{Mp)/y^^^'^ between the points in the static 
approximation and those at finite mass Mp to the mass of the B meson we obtain 
an estimate for fs- Here, however, there are two more sources of systematic errors 
which we include in our final result. The first is the variation on JeVMb in the 
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static approximation by using set 2 instead of set 1 for the central value. The second 
arises by choosing different sets of values for ^/a/Mp that are included in the 
interpolation. Our best estimate is obtained using set 1 and yJajMp = 0.21, 0.24 
and 0.26 as in the case of fu- Our result for /s is 

/b = 172 ± 24 (stat) li2(fits) li5(scale) MeV. (73) 

3/2 

We can now analyse the mass behaviour of $(Mp)/y^ in the continuum limit 
and deduce the leading correction in 1/Mp to the scaling law eq. (62). Figure 5 is 
the analogue of figure 1 at zero lattice spacing. If we parametrise <^{Mp) as 

$(M.) = $(Mp)^(l + ^ + ^) (74) 

we obtain C = 0.9 — 2.0 GeV. Thus, the conclusions reached at finite values of a are 
also valid in the continuum limit. The large uncertainty in C is partly due to the 

variation of ^(Mp)/ ^/a^^^ under changes in /Smin, and in particular to using cither 
set 1 or set 2 in the interpolation to the mass of the B meson (the highest values for 
C are obtained by using set 2 in the static approximation). 

The ratios foj Id and JeJIb are also of great interest. Since it can be ex- 
pected that systematic effects cancel partially in these ratios, one can attempt 
a straightforward extrapolation of the results quoted in refs. ^^'-^"^'■^-^^'-'^-^^'^-'^^ and 
refs. for /s^/Zs in the static approximation. In figure 6 we have 

plotted these extrapolations as a function of a^Ja, and one can sec that the results 
show a fairly flat behaviour as the continuum limit is approached. For our final 
results at a = we have again chosen /Jmin = 6.0 and continued the fits to lower 
values of /3 in figure 6 for illustrative purposes. We obtain 

^ = 1.08(8), ^ = 1.14(8) (75) 

JD JB 

^ = 1.07(5) t'l, (76) 

JB 

where we have used set 1 for the extrapolation of /b, / Jb in the static approximation. 
If instead one computes the ratio fosl fo using our estimates in eqs. (71), (72) one 
obtains foJfD = 1.08(14) '^sry good agreement with the above result, albeit 

with a larger error. 

3.5. Summary and comparison 

The analysis of lattice results for heavy-light decay constants shows that es- 
sentially a consistent picture emerges from all three methods of simulating heavy 
quarks. The errors are, however, still fairly large, so that this statement has to moni- 
tored as more precise lattice data become available. Adding all errors in quadrature, 
our best estimates for the decay constants fo, /d,, /b, /b, and Z^*^* are listed in 
tables 6 and 7. 
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Fig. 6. Extrapolations of fB^/fB and fos/fo to the continuum limit as a function of the string 
tension. The dotted lines are the continuations of the fits performed for /S > 6.0. Data at common 
values of a\/a have been combined in a weighted average. 



The large 1/M corrections to the scaling law for fp\fMp are now firmly es- 
tablished. In particular, the behaviour is not changed in the continuum limit. 
Furthermore, the heavy quark spin symmetry is manifest in the infinite mass limit, 
as U{M), the ratio of matrix elements of the axial and vector currents indeed ap- 
proaches unity. Also for this quantity the corrections in 1/M at the B and D meson 
masses are sizeable. 

The continuum limit requires further investigation. In particular it is necessary 
to increase the reliability of the extrapolations by going to smaller values of the 
lattice spacing, which correspond to (3 values larger than 6.3. Furthermore, as an 

important check against the inflticnce of lattice artefacts more data using improved 
actions and operators need to be computed in the range /3 > 6.0, in order to allow 
for a separate continuum extrapolation. So far, quite a number of results quoted at 
non-zero lattice spacing arc in agreement with the continuum result derived in this 
analysis (c.f. tables 1 and 2). In view of the large errors, however, this should not 
be taken as evidence that the continuum limit has essentially been reached already. 

We can now compare continuum lattice results for heavy-light deacy constants 
(also using the values quoted in other recent global analyses of lattice data is. 23, 25, 26 -j 
to the findings of other theoretical methods and experimental results. In tables 6 
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Table 6. Summajry of results for decay constants of charmed mesons. 



Source Author 



Lattice This work 



fo [MeV] fo^ [McV] 



Sum 
Rules 

Quark 
Models 



Exp. 



Martinelli , Flynn ^^'^6 
Sommer 

Neubert 

Dominguez 
Narisoii"4'"6 

Capstick+God£reyi5*5 
Oakesi57 

Hwang+Kim 
Hwang+Kim 

Average 



191 



+ 19 



-28 

205(15) 

164(20) 

170(30) 
139 - 234 
171(16) 

240(20) 



206 



+ 18 



-28 

235(15) 



174 - 269 
197(20) 

290(20) 



253(65) 265(68) 

< 310 241(21)(30) 
309(58)(33)(38) 



1.08(8) 
1.15(5) 



1.21(6) 
1.15(4) 

1.21(13) 
0.985 
1.045 
1.201 



and 7 we compare our best estimates with those from QCD sum rules, potential 
models and experimental results where available. The tables highlight the large 
uncertainties which are present in all determinations of these quantities so far. As 
mentioned before, lattice results suffer from systematic effects due to quenching, un- 
certainties in the extrapolation procedure and the treatment of heavy quarks. Esti- 
mates obtained using QCD sum rules are rather sensistive to the choice of input pa- 
rameters in the evaluation, which in the past has sparked a long controversy among 
different authors. More recent results from the sum rule approach 
however, agree quite well with those obtained on the lattice. It should be noted that 
QCD sum rule estimates for fpjfp are somewhat larger than the lattice values. It 
is, however, difficult to estimate the real uncertainties associated with the results 
obtained using QCD sum rules. An interesting result is quoted in ref. where 
an upper bound on fs, i.e. fs < 195 MeV is derived in a model-independent way, 
using rigorous QCD dispersion relations. 

Heavy-light decay constants arc hard to measure experimentally, since the lep- 
tonic width is small compared to the total width, and also because the presence 
of a neutrino in the final state complicates the analysis. Furthermore, the leptonic 
decay rate for a B meson is Cabibbo-suppressed, since it is proportional to IKjfcl^, 
which puts it beyond the reach of current experiments. For charmed mesons, only 
an upper bound can be quoted for /d ^'^^ (c.f. table 6), whereas leptonic decays of 
Dg are Cabibbo favoured, since their rate is proportional to Thus, fo^ has 

indeed been measured by several experiments 

The observation of large 1/M corrections to eq. (62) in lattice calculations has 
been confirmed by QCD sum rules Sum rule estimates for the 

coefficient C in eq. (74) vary in the range 0.7-1.3 GeV, in good agreement with 
the lattice result quoted in subsection 3.4. So far, only lattice calculations using 
NRQCD without the axial current normalisation have produced substantially larger 
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Table 7. Summajry of results for decay constants of bottom mesons. 



Source Author 



fB [McV] 



f§-' [MeV] fBjfB 



Lattice 



Sum 
Rules 



Quark 
Models 



QCD DR 



+27 



This work 172 
Martinelh^s, Flynni^'^s 
Sommer ^'^ 

Neubert 

Elctskii+Shuryak 

Bagan et al. '^^'^ 

Dominguez ^^'^ 
Narisoni44'"6 

Capstick+Godfreyi^^ 
Oakcs 1" 
Hwang+Kim 
Hwang+Kim 

Bovd (-/ al. < 195 



245 



+50 
-47 



-31 

175(25) 

180(46) 276(37) 

190(50) 200 - 300 
150- 175 165- 200 
195 - 245 

133 - 183 

209(34) 259(41) 
155(15) 

201(51) 



1.14(8) 
1.15(5) 
1.10-1.18 



1.22(2) 
1.16(4) 

1.35(18) 
0.989 
1.053 
1.173 



corrections, such as C = 2.8(5) GcV quoted in rcf. 

To summarise, we have shown that lattice calculations using unimproved conven- 
tional and static heavy quarks have produced results for heavy-light decay constants 
with a total accuracy of 15% for charmed and 30% for bottom mesons. A number 
of systematic errors necessitate further studies, as signified by the rather large er- 
rors. Preliminary results from runs with dynamical quarks ii^.iae g^ggest that decay 
constants could increase by up to 10% in the unquenched theory. However, unlike 
some time ago, there is now agreement between lattice calculations and QCD sum 
rules, and also the few experimental results. It is of great importance to increase 
the precision of lattice data, especially in view of future dedicated experiments for 
B physics. 



4. S°-BO Mixing 

We will now turn the discussion to lattice calculations of the B parameter Bb, 
which is relevant for B^-B^ mixing. The lattice estimates for Bb and the decay 
constant /b from the previous section can then be combined into the quantity 
fss/BB, which is of particular interest, since it constrains the possible values of 
the parameters p and -q in the standard Wolfcnstein parametrisation of the CKM 
matrix. At the end of this section, we will therefore attempt to quote a common 
estimate for this quantity using data from several lattice calculations and study its 
phenomenological implications in section 5. 



4.1. The AB = 2 four-fermion operator on the lattice 



In the continuum the interpolating operator for oscillations between a i?° and a 
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meson is the four-fermion operator Ol defined by 

Ol = (67^(1 - 75)5) (6741 - 75)?). (77) 

This operator enters the AS = 2 effective Hamiltonian . The ampUtude of B^-B^ 
mixing is usually expressed in terms of the B parameter, which is the ratio of the 
operator matrix element to its value in the so-called vacuum insertion approximation 

o _ (B^\Ol{^^)\B'^) _ {W\OUl^)\B'^) 

Here, ^ is a renormalisation scale, and is the temporal component of the heavy- 
light axial current Af^ = bj^j^q. The dependence of Bb on the renormalisation scale 
can be removed by multiplication with a factor derived from the anomalous dimen- 
sion of Ol. At leading order (LO) or ncxt-to- leading order (NLO) respectively, a 
renormalisation group invariant B parameter can be defined by 



B}f = a,{n)-^'^-BB{n) (79) 
= as{,.)-'l^^{l + ^Jr,,)Bs{„), (80) 

where /3o = 11 — 2n//3 and J„j, is derived from the one- and two-loop anomalous 
dimensions 7}°^ and 7^^^ of the operator Ol- In the MS scheme one obtains 



7?' =4, 7?' = -7+ in,. (82) 

With these definitions is renormalisation group invariant only up to leading 
order and next-to-leading order, respectively. 

If the operator matrix element is evaluated on the lattice using Wilson fermions, 
the contributions of the corresponding operators with opposite and mixed chirality 
need also to be taken into account, as was explained in subsection 2.3. Since every 
operator has to be renormalised separately, the matching between the lattice oper- 
ators and their continuum counterparts is potentially an intricate problem. Hence, 
in the case of four-fermion operators, the importance of a non-perturbative deter- 
mination of the various matching factors is even greater. Attempts in this direction 
have been reported in the context of the B parameter Bk 

Lattice results for the B parameter Bb have been published recently using the 
static approximation ^"^^'^^^'^^^ and also propagating heavy quarks Ear- 
lier attempts are discussed in refs. I09'i26,i35 

Here wc shall describe in some detail the evaluation oi Bb in the static approx- 
imation and the subsequent matching of the results in the effective lattice theory 
to the AS = 2 Hamiltonian in the continuum, which plays the role of the "full" 
theory, even though it is an effective theory in itself. The matching is usually per- 
formed in two steps (although there is in principle no need to do so), in which one 
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first matches the operator in full QCD at scale fj, = nib to the continuum effective 
theory (i.e. HQET). In the second step one matches the basis of local operators 
in the continuum effective theory to those in the effective theory on the lattice. 
The following diagram illustrates the two-step process and lists the operators and 
renormalisation scales at every stage: 



Lattice 

static approx. 
scale 

Ol, Or, On, Os 



Continuum 




Continuum 




AS = 2 effective 


step 1 


HQET 


step 2 


Hamiltonian, 




scale fi 




OL{mb) 









Thus, when relating the continuum full theory to the continuum effective theory 
the operator Os is generated at order ag, owing to the mass of the heavy quark 

Os={b{l-^,)q){b(l-7^)q). (83) 

Furthermore, due to explicit chiral symmetry breaking, the following operators mix 
with Ol in the discretised theory 

Or = (%(l + 75)<z) (67^(1 + 75)9), (84) 
On = (67^(1 - 75)g) (67m(1 + 75)9) , 

+(^,.(1 + 75)9) (h^i - 75)9), 

+2(6(1-75)9) (&(1+75)9), 

+2(6(1 + 75)9) (&(1- 75)9)- (85) 

The relation between the B parameter in the full theory and the relevant operator 
matrix elements in the effective lattice theory is then given by 

ExZr'lB^ Ox{a) B^) 

BB(mb) = ^ 5 — X = L,R,N,S (86) 

l{zr'{0\A4B0)f 

Here the matching factors Zjj** and Zj** for the axial current are the products of 

the renormalisation factors required for steps 1 and 2, which arc chosen such that 
all dependence on the scales /U and a in Bsimb) vanishes. So far, only perturbative 
estimates exist for these matching factors, and we refer the reader to appendix A 
for explicit expressions. Here, we only mention that a complete matching at order 
as can be performed, since the two-loop anomalous dimensions for step 1 have been 
computed in the MS scheme 

4.2. Lattice results for Bb and Bb^ 

Results for the B parameter Bb are compiled in table 8. In order to facilitate the 
comparison of different calculations we have chosen to quote BB{mb) at a reference 
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Table 8. Results for the B parameter Bs(mi,) obtained in the quenched apijroxiriiation. B^*^ and 
B^^'-' are the renormalisation group invariant B parameters obtained from B g (mj ) as explained 
in the text. Due to our particular choice of parameters, the results for Bg may differ slightly 

from those quoted in the original articles. The static approximation was used in refs. 100,162,16,'!^ 
Furthermore, UKQCD as well as the authors of ^^'^ use an improved action (csw = 1) for the light 
quark sector. 



Collab. 


P 




Bb {m,,) 




jjNLO 




6.0 


0.98(4) 

M = 4.33 GeV 


0.97(4) 


1.42(6) 


1.54(6) 


UKQCD 


6.2 




0.69 +1 tl 


1.02 tl tl 


1.10 ^2 


G+M162 


6.0 




0.63(4) 
0.73(4) 


0.92(6) 
1.07(6) 


1.00(6) 

1.16(6) 


JLQCDi" 


6.3 
6.1 




0.840(60) 
0.895(47) 


1.23(9) 
1.31(7) 


1.34(10) 
1.42(8) 




"oo" 


0.96(6)(4) 
^ = 2 GcV 


0.89(6)(4) 


1.30(9)(6) 


1.42(10)(6) 


ELC 


6.4 


0.86(5) 

M = 3.7 GeV 


0.84(5) 


1.24(7) 


1.34(8) 


BDHS 


6.1 


1.01(15) 
M = 2 GeV 


0.93(14) 


1.36(20) 


1.48(22) 



scale TOft = 5 GcV, although some authors prefer a different value. For fj, < 5 GeV, 
this conversion has been performed at leading order according to 

BB{mt = 5 GeV) = Bsif^) (^4^) ■ (87) 

We have used the expression for as{\£) at leading order, viz. 

47r 

= /?oln(MVA-) ' ^''^ 
which was evaluated for A = 200 MeV and n/ = 4 active quark flavours. Requiring 
as to be continuous at the threshold /x = rUb, we find a^°(m{,) = 0.2342. The value 
of the renormalisation group invariant B parameter B^ in table 8 was obtained 
from eq. (79) using the value of a^°(m6). Similarly, B^^ was obtained from 
eq. (80) by inserting the value of (xf^'^{ra],) = 0.1842. 

By comparing the results for BB{mb) in the table one observes a cluster of values 
in the range 0.80 — 0.90, as well as a few results clearly above and below, despite 
the fact that similar lattice spacings have been used. As was first noted in ref. 1°°, 
ambiguities in the perturbative matching procedure are largely responsible for this. 
As the matching factors Z^^*', Z^^^, X = L,R,N,S in eq. (86) are themselves 
products of matching factors known to 0{as), the ratios Z^'^^ / {Z^'^^)'^ contain some 
contributions of order a^. By including or excluding these next-to-next-to-leading 
order terms, the following three possibilities can be applied: 
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Table 9. Results for BB{mt,) in the static approximation for matching procedures (Mi) — (M3), 
using results for the individual matrix elements from different groups. The couplings as(^), as (mi,) 
were evaluated using rnf, = 5GcV, nj = 4, A^^^ = 200 McV. Only statistical errors arc shown. 
Tadpole improvement has not been applied, except that a mean field improved coupling ay (5*) 
at aq* = 2.18 has been used. 



method 


as{mb) 


UKQCD 


G+M ^'^^ 


Ken '"'^ 


(Ml) 


0.2342 


0.81(4) 


0.82(3) 


0.87(3) 




0.1842 


0.80(5) 


0.81(3) 


0.85(3) 


(M2) 


0.2342 


0.60(4) 


0.59(3) 


0.42(3) 




0.1842 


0.64(5) 


0.64(3) 


0.51(3) 


(M3) 


0.2342 


0.87(5) 


0.83(2) 


0.96(3) 




0.1842 


0.88(5) 


0.85(3) 


0.98(3) 



(Ml) include O(a^) contributions by using the ratios ZJ^^'Vl^yi''*)^- 

(M2) exclude 0{al) contributions in the individual Z- factors by expanding Z^^* 
and Z^^* separately to order a^. 

(M3) exclude 0{al) contributions by expanding the ratios Z^'^^ / {Z^'^^)'^ to or- 
der ag. 

In ref. ""^^^ it was suggested that the low values for BB{mi,) reported by UKQCD 
were due to using procedure (M2), and that therefore the large discrepancy between 
refs. and could largely be explained by the amplification of ambiguities by 
employing a particular matching procedure, rather than the different values of Csw 
or (3 in the two simulations. In view of this issue, we wish to test the consistency 
of results from different groups, by listing in table 9 the values for BB{mb) from 
refs. 100,162,163^ which were obtained by applying methods (Mi), (M2) and (M3) for 
uniform choices of input parameters. Indeed it turns out that method (M2) system- 
atically produces lower values of Bsirni,). However, for each common procedure 
the discrepancy between the results for Cgw = (ref. ^^'^) and Cgw = 1 (refs. ioo>i62-j 
cannot fully be resolved, so that other effects must also contribute to the observed 
differences. 

One concludes that, although different prescriptions should be equivalent in 
principle, there are discrepancies of up to 25% in the final values of BB{mb) for 
Csw = 1 and even larger variations for unimproved Wilson fermions. As each of 
the methods (Mi) — (M3) is perfectly justified, these differences may be viewed as a 
systematic error on Bb as a result of perturbative matching. As far as lattice results 
for Bb in the static approximation are concerned, these ambiguities do not arise 
exclusively from the discretisation of the theory, since the matching step between the 
continuum full and continuum effective theories contains a large 0{as) contribution 
(see also appendix A). This problem can only be adressed in a two- loop calculation 
for this matching step. 

Different procedures should eventually yield the same answer in the continuum 
limit, provided that lattice artefacts are not so large as to spoil the extrapolation 
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a 0. Clearly, the currently available data for Bb do not allow for such an 
extrapolation as in the case of decay constants discussed in the previous section. 
However, preliminary results in the continuum limit for Bb and some SU(3)-flavour 
breaking ratios involving the B parameter have been reported Also, by 

comparing the results from rcfs. i"0'ifi2^ which were both obtained for Cgw = 1) one 
observes only a weak dependence on the lattice spacing a. 

In the absence of further data which might clarify the issues of matching as well 
as the influence of lattice artefacts and the size of 1/M corrections to the results in 
the static approximation, we resort to quoting a central result for Bb with an error 
which encompasses the observed variations. Our best estimate is 

BB{mb)= 0.85 til, (89) 

which corresponds to 

B]f = 1.24 111 B^^° = 1.35 til, (90) 

so that our final estimate for the renormalisation group invariant B parameter can 
be summarised as 

Bb = 1.3 tl (91) 

These numbers can be compared with the results from QCD sum rules. The authors 

of ref. "'^^^ obtain BB{5GeV) = 0.99(15), which, in view of the large errors, is in 
agreement with our estimate. An earlier attempt to calculate the matrix element 
of Ol in the sum rule approach resulted in a very low value 

We can combine our value for Bb with the result for the decay constant /s in 
the previous section. Using /s = 172 tg^MeV from table?, we find 

/bV^=195 t4oMeV. (92) 

This is one of the central results of this review and can be taken as the present 
estimate for this quantity treating as many individual results as possible on an 
equal footing. 

Now we turn the discussion to SU(3)-flavour breaking ratios involving the B 

parameter and the decay constant. As becomes clear from table 10 the B parameter 
shows very little dependence on the mass of the light quark, quite in contrast to 
the decay constant. The results for BbJBb^ can therefore be summarised by the 
value 

BbJBb, = 1.00(2). (93) 

The result can be combined with our estimate of fBjfB^ = 1.14(8) in table? and 
eq. (?5) to compute the ratio 



^sd = = 1-30(18). (94) 
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Table 10. Results for the SU(3)-flavour breaking ratio Bg^/Bg^ and the ratio defined in 

eq. (94). Note that Bg^/Bg^ = Bg^{nit,)/Bg^(mi;), and thus this ratio is independent of the 
details of the conversion to the renormalisation group invariant B parameter. 



Collab. 


P 






Comments 




Kon 


6.0 


0.99 ti(l) 









G+M 


6.0 


1.01(1) 


1.38(7) 


S t tit IC ^ w — 


1 


UKQCD 


6.2 


1.02(2) 


1.34 tl 11 


StcltlC^ Csw — 


1 


JLQCD114 


6.3 


1.05(3) 




conv., Csw = 







6.1 


0.99(3) 




conv., Csw = 





BBS 


"oo" 




1.49(13)(31) 


conv., Csw = 





this work 






1.30(18) 


eqs. (75,93) 





More results for which have not been used in the derivation of this result, 
are listed in table 10. It is worth noting that apart from the estimates for in 
the static approximation (ref. i00'i62)^ there is also a preliminary result lo^ which 
treats heavy quarks in the conventional approach. Moreover, in that study was 
computed directly from the ratio of operator matrix elements, viz. 

_ (B;|Oi|b;) Ms. 

over a range of lattice spacings. The results in table 10 and their statistical errors 
can be combined in a weighted average to yield 

= 1.38(15), (96) 

where we have added in quadrature a systematic error estimated from the spread 
of central values in the table. 



5. Analysis of CP violation 

In this section we shall combine lattice results with experimental data and other 
theoretical estimates to place constraints on the most poorly known elements of the 
CKM matrix and on the unitarity triangle. 

5.1. Parametrisations of the CKM matrix and the unitarity triangle 

We start the discussion by recalling the basic definitions and relations involving 
the CKM matrix. As is well known, tlw^ CKM matrix Vckm relates the gauge cigc^n- 
states appearing in the SM Lagrangian to mass eigenstates. For flavour-changing 
charged current transitions this has the consequence that, in addition to the domi- 
nant transitions between up- and down-type quarks, u ^ d, c ^ s and t ^ b, there 
are further transitions of lesser strength. The corresponding Lagrangian describing 
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flavour- changing charged current interactions in the hadronic sector has the form 

C^^S = - ^(^' ^' - VcKM I s \Wl+ h.c., (97) 
where the CKM matrix 

VcKM = I Vcd Vcs Vcb I (98) 




Vud 


Vus 


Vub 


Vcd 


Vcs 


Vcb 


Vtd 


Vts 


Vtb 



is a unitary matrix in flavour space. Since the dominant transitions are proportional 

to the diagonal matrix elements Vud, Vcs and Vtb, one expects a hierarchical structure 
of VcKM- An approximate parametrisation of the CKM matrix which takes account 
of this structure is due to Wolfenstein^^^. By expanding Vckm in powers of the 
Cabibbo angle |V^s| = A ~ 0.22 to order A^ one obtains 

/ 1- f A AX^p-i7^) \ 

Vckm ^[ -A 1 - f AX^ , (99) 

\ AX^{1- p-irj) -AX"^ 1 / 

with the remaining parameters A, p and 1] of order one. 

In the case of three generations of quarks and leptons, Vckm can be parametrised 
in terms of three angles and one observable complex phase, of which the latter is 
required to describe CP violation. As has been pointed out in ref. a consistent 
treatment of CP- violating effects necessitates the inclusion of higher terms in the 
Wolfenstein parametrisation, viz. 



Vckm ^ [ -A - iA'X^r] 1 - ^ AX^ \, (100) 



1-^ A AX^p-iTj) 

AX^{l-p- irj) -AX^ - iAX^T] 1 
where the rescaled parameters p and rj are expanded as 

-p = p{l-^ + 0(A^)), ^ = ^(1 - ^ + 0{X')). (101) 

The relative size of CKM matrix elements is easily recognised in the Wolfenstein 

parametrisation: the diagonal elements are of order one, and \Vus\, \Vcd\ are both 
of order 20%. The relative size of \Vcb\ and \Vts\ are 4%, whereas \Vub\ and \Vtd\ 
are of order 1 %. 

The Wolfenstein parameters A and A are rather well determined experimen- 
tally "9.i70,i74,i76 

A = IKsl = 0.2205 ±0.0018 (102) 
A = =0.81 ±0.058, (103) 



y2 

^ US 
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whereas the elements \ Vub\ and \Vtd\ have uncertainties of around 30%, so that the 
parameters p and ry are rather poorly known. 

The unitarity of Vckm imposes the following conditions on its elements 

Vij V*, =0, k. (104) 

There are six such conditions, each of which can be represented graphically by a 
triangle. The triangle relation for the most poorly known CKM matrix elements 
reads 

Vudv:, + KrfV;; + Vt^v;, = o, (io5) 

and its representation is shown in figure 7. 

In order to constrain the values of p and rj one needs information on \ Vtd\, which 
can be extracted from B^-B^ mixing. The expression for the mass difference AMd 
between the B^ and B^ states reads 

AM, = -^-^r,B,S{mt/Mw)fB,BB, \VM^ (106) 

where Gi? is the Fermi constant, and Myy is the mass of the VF-boson. The quantity 
TjBd = 0.55±0.01^^^ parametrises QCD short-distance corrections, and S{mt/Mw) 
is a slowly varying function of rrit/Mw In the following we shall evaluate 

the above expression using our lattice estimate for fB^\J Bb^^- 

5.2. Constraints on the Wolfenstein parameters p and r] 

Using lattice data in conjunction with experimental data and other theoretical 
estimates, we now want to analyse B^^-B^ and B'^-B^ mixing in order to place 
constraints on the parameters p and rj discussed in the last subsection. 

As is obvious from eq. (106), the quantity fs^ \J Bb^ plays a central role in this 
study. Further input parameters are the running mass of the top quark, TO((to(), the 
short-distance corrections rjB^ and the experimentally measured value of AM,. For 
B^- B^ mixing, the quantity is also required. A compilation of input parameters 
for this study is listed in table 11. 
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Solving for \Vtd\ in eq. (106) and inserting the relevant input parameters as well 
as r]B = 0.55±,0.01, one finds 

|V^d| = (8.94^^i|)xlO-^ (107) 

and thus the relative uncertainty in \Vtd\ according to this estimate amounts to 
15 — 20 %. The results for \Vtd\ can be translated into a constraint on p and rj via 
the unitarity triangle 

where we have used the present determination of \Vcb\ 

\Vcb\ = 0.0393±, 0.0028. (109) 

Further constraints on p and rj can be obtained from indirect CP violation in the 
K^-K^ system through a theoretical analysis of the parameter ck ^''^'^^'^'^^i. This 
requires as additional input an estimate for the kaon B parameter B^- Several 
theoretical methods, including lattice simulations/^^' quote values for Bk 
which are compatible with the range 

= 0.75 ±0.10. (110) 

There has been considerable progress in calculating Bk on the lattice However, 
a number of systematic effects such as the continuum limit, the chiral behaviour of 
Bk computed using Wilson fermions, and the dependence on the number of flavour 
require further studies A more detailed discussion of lattice results for Bk 

is, however, beyond the scope of this review. 

Herrlich and Nierste have obtained constraints on p and rj from a theoretical 
analysis of ck using the complete AS' = 2 effective hamiltonian at next-to-leading 

order. Using very similar estimates for /s^ Bb^ and (i.e. /s^ \f^Bd = 200 ± 
40MeV, .^^^ = 1.32 ± 0.12) compared to this review, they find 

-0.20 < p < 0.22, 0.25 <^< 0.43, (111) 



Vtd 



cb 



1.03 



+0.18 
-0.23' 



(108) 



Table 11. Input parameters for the study of p and rj. 
Quantity Value Source 

fsa \Jbb^ 195 +4oMeV Lattice, this work 

1.38 ±0.15 Lattice, this work 

Mb^ 5279 MeV Exp. 
Mb, 5369 MeV Exp. 
AMd 0.464 ± 0.018 ps'^ Exp. 
mt(mt) 165±9GeV Exp. ^^^'^^^ 
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in good agreement with a similar study by Ali and London 

Now we turn our attention to B^-B^ mixing. The ratio of mass differences for 
the ByB^ and B^-B^ systems is 

AM, _ Pb^Bb^ Mb, [y^.P , . 

AM, Vb,P^Bb,Mb,\V,,\^- ^ ' 

Since riBs = Vsa this becomes 

AM, ^ Mb, \Vts\^ 



AMd ^'"Mb.IVuI^ 



= (1.40±0.15)|U, (113) 

where we have used the parameters Usted in table 11. This result can be used 
together with the measured value of AMd to predict AM,, provided that the allowed 
range of |T4sP/|Vtdp is determined. This has been performed by Ali and London 



who have analysed |Vtc;|/|Vts| as a function of /b^ \J Bb^, assuming that B^ = 



0.75 ± 0.10. Using our estimate of /s^ \J Bb^ = 195 llo^eV one reads off the 
allowed range of |Vtd|/|V(s| from figure 4 in ref. ^"^^ 



\Vtd\ 



0.228 ± 0.040 (114) 



\Vts\ 

where the first error is duo to the uncertainty in Bk, and the second reflects the 



error on Jb^ y Bb^ ■ This yields 

AMs = 12.5 ± 0.48 ± 1.3 ± 4.4 t|:^ps"^ (115) 

Here the first error is due to the experimental error on AM^, the second due to the 
error in and the third and fourth to the uncertainties in the matrix elements 

Bk and /s^ y^B^^, respectively. Adding the errors in quadrature one gets 

AMs = 12.5t^:oPS-^ (116) 

Using the updated results for the Bg lifetime, tb^ = 1-61 ± 0.10 ps^^^, this corre- 
sponds to a value for the B^-B^ mixing parameter Xg of 

Xs = 20.1^^1:^. (117) 

Our estimates for AMs and Xs can be compared to the experimental lower bound 
of 

AMs>7.8ps-\ a;^>12.6 (95% C.L.) (118) 

quoted by the ALEPH Collaboration This indicates that, although the 
experimental lower bound is still smaller than the central value of the SM prediction, 
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experiments start to exclude parts of the allowed range of AM^. Therefore, exper- 
imental data for AM4 and AMg can further restrict the allowed range of (p, rj), a 
fact that has already been exploited in deriving the constraints in eq. (Ill) 

The constraints on CKM parameters have sharpened significantly in recent years, 
not least thanks to more accurate and consistent determinations of weak matrix 
elements using lattice techniques. 

6. Concluding Remarks 

In this review we have combined a variety of lattice results for heavy-light decay 
constants and the B parameter Bb into a common estimate for /bVBb- A central 
part of this study was the analysis of the continuum limit of decay constants and the 
assessment of systematic errors. Our result for JbVBb has been combined with 
other theoretical and experimental input to constrain some of the CKM matrix 
elements. Although the overall errors on lattice results in the heavy quark sector 
are still rather large, the progress in lattice calculations over the past years has 
helped a lot in order to sharpen the bounds on CKM matrix elements and the 
unitarity triangle. 

The interplay between different theoretical tools plays an important role in this 
analysis, as it serves to extract a consistent picture in situations where each method 
has its own limitations. An example for this is the study of 1/M corrections to the 
scaling law for decay constants. Lattice simulations arc currently limited by a 
number of systematic effects, such as quenching, which cannot fully be addressed 
at present due to lack of resources. However, in order to increase the precision of 
lattice calculations, not only larger computing power but also new theoretical ideas 
such as non-perturbative improvement, arc of great importance. 

Here we have concentrated on leptonic decays of heavy quarks. This is by 
no means the only area where lattice calculations have made an impact on the 
determination of CKM matrix elements. Indeed, another poorly known element, 
Vub, can be determined by considering semi- leptonic B decays, such as _B ^ irli/i or 
B — » pivi- Lattice simulations offer model-independent estimates for the relevant 
form factors, and a variety of calculations has already been performed Likewise, 
lattice results for rare decays such as _B ^ K*"/ have been reported in the past 
few years. All of these calculations serve ultimately to test the consistency of the 
Standard Model, and it is evident that the lattice is a versatile tool in these analyses. 

With the advent of dedicated experimental facilities for B physics, such as 
HERA-B, BaBar and LHC-B, together with an improved theoretical understand- 
ing of heavy quark physics, the scene is set for an exciting new period in particle 
physics. 
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Appendix A. Matching Factors for the Static Approximation 

In this appendix we list the expressions for the perturbative matching factors 
for the heavy-light ax;ial current and the AB = 2 four-fermion operator in the static 

approximation. As described in subsection 4.1, the matching is usually performed 
in a two-step process, where one first matches the full theory in the continuum to 
the effective continuum theory, which is then matched to the effective theory on the 
lattice. 

Throughout this appendix we use the following convention for the one- and 
two-loop coefficients of the perturbative /3 function 

2 38 
/3o = 11 - ^Uf, /3i = 102 - —Uf, (A.l) 

where n/ is the number of active quark flavours. 

Appendix A.l. The axial current 

A full 0{as) matching between the axial current in full QCD and its coun- 
terpart Afj, in the continuum effective theory can be performed, since the two-loop 
anomalous dimension of has been computed For typical hadronic scales 
fj, < nib, and using the naive dimensional regularisation scheme one obtains ^'^'^^ 



where 
and 



C2 = -8/3 (A.3) 




(0) 

^^1,^ = 1^' Ja = ^{^-^]. (A.4) 

7f = -4, 7^' = - + f (^-5) 

The general expression for the matching of the current to the current A^ in the 
eff'ective theory on the lattice is 

= (1 + [2 ln(M'a2) + Da] ) A^. (A.6) 
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Numerical values of Da depend on the choice of Csw and are tabulated in table A. 1. 

Eq. (A. 6) is usually evaluated using a tadpole improved or boosted value for 
the lattice coupling a|f**, for instance 5s = a^f^/^O' where uq is the average link, 
uo = (|ReTrP)^/^, or alternatively uq = l/(8Kcrit)- Combining measured values 
for Uq with its perturbation expansion, viz. 

uo = l + u'^o^as + O(a^) (A.7) 

one can write down the tadpole improved version of eq. (A. 6) 

^;. = vW(l + ^[21n(MV)+2)^-i«W]) V (A.8) 

The square root of uq has to be taken since the static-light axial current contains 
only one relativistic quark field. The matching factor between and could in 
principle also be determined non-perturbatively using the method in ref. . Com- 
bining the factors for the two matching steps one finds 



-'A 



as{mb) \ '"^ f asif^) - as{mb) % as{mb) \ 

asili) J \ 47r ^ 3 47r j 

X (l + ^ [2\n{n'a') +Da- i^'^]) • (A.9) 



Appendix A. 2. The four-fermion operator 

The perturbative matching for the four-fermion operator is more complicated 
due to operator mixing. Apart from the operator Ol, the basis of local operaotrs in 
the continuum effective theory consists also of the operator Os defined in eq. (83). 
The full 0{as) matching relation between and the operaotrs and Os in the 
continuum effective theory is given by 




as{mb)Y'\ 72? asjmb) ^ ~ , , 

722 111 

^Y^B^Osi,), (A.10) 
as(/i) / 47r 

where in naive dimensional regularisation^^ 

Bi = -14, B2 = -8 (A.ll) 

and 64,66-68 



Leptonic decays of heavy quarks on the lattice 51 



Table A.l. Values for the coefficients Da, Di^, Dn and Dat for unimproved and improved Wilson 
fermions according to refs. In the case of the four-fermion operator for Csw = 1, the 

conventions of ref. "^"^ are used. 





Da 


Dl 


Dr 


Dn 





-27.17 


-38.91 


-1.61 


-14.40 


1 


-20.20 


-22.50 


-5.40 


-14.00 



' - \ 7^? 7^? J " U/3 -8/3 ) ' (^-^3) 

(1) 808 52 2 64 
7li^=-^-^-^ + yn,. (A.14) 

When matching Ol(m) to its lattice counterpart Ol(p), one has to consider the 
lattice operators Or{^) and On{ii) (eqs. (84,85)), such that 

Oi(M) = (^l + ^^^[41n(/iV)+£>i]^ Oi(a) 

+^£^i5^0^(a) + ^£^2?^0^(a). (A.15) 
47r 47r 

Since Og is generated at order one finds 

Osin) = ds{a) + higher orders. (A. 16) 

Values for the coefficients Dl, Dr and D^v have been determined for Cs„ = 
and 1 64,62,77 g^j.g iig^;gj]^ in table A.l. We have chosen the convention that the 
static-light meson propagator is proportional to e~^*, such that the reduced value 
of the quark self-energy has been used^^ in the evaluation of Da and D^. 

Combining the two matching steps and using tadpole improvement, we find the 
following expressions for the renormalisation factors Zjj^*, X — L, R, N , S 



Oisip) J y 47r 47r 

xt^o (l + ^ [41n(M^a2) + - 



(-2)^f^ (A.17) 



^ as{p) J y ' 47r " " A-K 

^^[Un{^Jt''a')+Dn] (A.18) 
47r 

_^stat ^ j a«(mb) y^ ^ a^(/i) - as(m6) ^ ^^a^(mb)' 



— T-r \ ^ 'J — ■'■4 
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X 



47r 



(A.19) 



'Stat 
S 




4tt 



(A.20) 



One can now use the above expressions directly in the calculation of BB{mb) or, 
alternatively, perform a systematic expansion of Z^*** and Z^^^, X = L, R, N, S 
to order a^. Since the B parameter is a ratio of matrix elements, another possibility 
is an expansion of the ratios Zjj'^'/l-^'^''')^, X = L, R, N, S to order ag. 

Appendix B. Continuum Limit in the Light Quark Sector 

Here we describe the continuum extrapolation of quantities in the light quark 
sector, which are commonly used to set the scale in lattice gauge theories. We 
will focus on the rho meson mass Mp, the pion decay constant and the IP-IS 
splitting in charmonium, Aip_is. 

For these three quantities the question arises whether or not they can be com- 
puted with high enough precision and reliably be extrapolated to the continuum 
limit: 

• precise lattice estimates for require knowledge of the current normalisation 
constant Za, which for most of the available lattice data is only known in 
perturbation theory. 

• Mp is usually obtained with small errors, but studies with high statistics 
show oscillations in the relevant correlation functions, which, despite the fact 
that they can be modelled, makes it more difficult to extract a precise mass 
value. Furthermore, as has been discussed in ref. the a dependence of 
Mp computed with unimproved Wilson fermions is very strong, such that the 
continuum extrapolation may not be too reliable. 

• the splitting Aip_is shows only a weak a dependence. However, the estimates 
published so far have still relatively large statistical errors 

This discussion underlines the importance of obtaining low energy quantities with 
higher precision, and this may now be possible using the recently proposed non- 
perturbatively improved Wilson action which is expected to yield lattice data 
that are completely free of lattice artefacts at order a. Since also the current 
normalisation constants have been determined non-perturbatively , one can expect 
highly accurate continuum results for as well. For the time being, however, 
we must restrict our attention to the unimproved Wilson lattice data for 

improved actions arc not yet available over a large enough range in /3 to allow for a 
continuum extrapolation. 

The extrapolation of hadronic quantities in the light quark sector proceeds by 
forming dimensionless ratios, say /tt/Vo"; using lattice data for the string tension, 
and extrapolating it to the continuum limit as a function of ay/a. In table B.l 
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Table B.l. Lattice results for the string tension a^/a. 



a 




Author 


6.8 


0.0730(12) 


BaU&SchiUing^y'' 


6.5 


0.1068(10) 


UKQCD 


6.4 


0.1215(12) 


BaH& Schilling 




0.1218(5)(23) 


UKQCD 




0.1216(11) 


combined 


6.2 


0.1619(19) 


BaU&SchiUingi96 




0.1609(28) 


UKQCD 




0.1608(7)(16) 


UKQCD 198 




0.1612(12) 


combined 


6.0 


0.2265(55) 


Bah&SchiUingi96 




0.2182(21) 


Michael & Perantonis 




0.2154(10)(40) 


UKQCD 198 




0.2185(18) 


combined 


5.9 


0.2702(37) 


MTc20i 


5.8 


0.3302(30) 


MTc^oi 


5.7 


0.4099(24) 


MTc^oi 



we list the lattice data for the string tension that have been used throughout this 
review. In cases when several estimates are available at one /? value they have been 

combined in a weighted average. Estimates of a^/a at "intermediate" values of /3 
(e.g. (3 = 6.17) were then obtained by linear interpolation of ln(a-^/a) as a function 
of (3. 

The conceptual advantages of the hadronic scale ro described in ref. make it an 
ideal quantity to perform the extrapolation to the continuum limit. A compilation 
of lattice results for ro/a can be found in table B. 2. Since ro/a is so far only available 
for (3 > 6.0, the Wilson data at lower (3 values (in both the light and heavy quark 
sectors) cannot be used in the extrapolations to a = 0. Thus, we have decided 



Table B.2. Lattice results for the hadronic scale ro/a. 



p 


ro/a 


Author 




6.5 


11.23(21) 


UKQCD 195,197 




6.4 


9.90(54) 


Bali & Schilling 


196,197 




9.75(17) 


UKQCD 198 




6.2 


7.38(25) 


Bah & Schilling 


196,197 




7.29(17) 


UKQCD 198 






7.27(3) 


SESAM202 




6.0 


5.44(26) 


Bali & Schilhng 


196,197 




5.47(11) 


UKQCD 198 






5.35 tg 


SESAM202 
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to use the string tension, keeping in mind that eventually ro should be used when 
more precision results for hadronic quantities become available. 

In order to obtain /tt/v^ and Mpj^ in the contimmm limit, wc shall use 
the results quoted by the GFll collaboration Their simulations were in- 
tended as a comprehensive study of the continuum and infinite volume limits of the 
light hadronic sector in the unimproved Wilson theory. There are, in fact, many 
more lattice results for and M/8'99'i02,i09,ii6,i28,i33,207-2ii^ ^jg^ improved 

actions and mostly with small statistical errors. However, a reliable 

continuum extrapolation of all available data is diflacult to perform without a con- 
sistent treatment of systematic errors. This is the main reason why wc have chosen 
to use only GFll's results in the present study. The accuracy and reliability of the 
extrapolated values can be checked by comparing to other simulations. 

In figure B.l (a) we show the continuum extrapolation of /tt/v^- For the axial 
current normalisation factor Za we used the tadpole improved pcrturbative expres- 
sion evaluated using the boosted coupling g\, — Qq/uq. It is seen that most other 
simulations are in fact consistent with GFll's data. In principle the error on the 
extrapolated result could be reduced by including more data points. However, this 
would imply a very selective use of the available data, since overall compatibility of 
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0.285(13). 



(B.l) 



We have checked the stabihty of this result by using different prescriptions to evalu- 
ate Za and find that the resulting variation is much smaller than the quoted error. 

Figure B.l (b) shows the corresponding extrapolation of Mpj^fa. The result at 
a = is 



The figure illustrates the strong a dependence of this ratio in the unimproved Wil- 
son theory (note that Mpj %Ju is plotted on a much smaller scale than /tt/ \fo). It is, 
however, reassuring to notice that more recent simulations using a non-perturbative 
value of Csw obtain value's for Mp/^/a which arc consistent with the extrap- 

olated Wilson result, whilst showing very little residual a dependence. 

In order to compute Aip_is/v^ continuum limit we have extrapolated 

the results quoted in ref. and obtain 



It should be emphasised that the non-pcrturbatively improved action could be 
used to obatin much more precise estimates of the above quantities. One will thus 
be able to set the scale much more reliably in future simulations. 
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